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INTRODUCTION 


With  the  increasing  cost  of  asphalt  over  the  past  decade,  it  has  become 
more  economically  attractive  for  agencies  such  as  the  U.  S.  Air  Force  to  recycle 
existing  airfield  pavements  instead  of  overlaying  them  with  new  asphaltic  con¬ 
crete  pavement.  Several  recycling  methods  are  available,  including  some  that 
apply  heat  to  the  surface  of  the  pavement  thereby  increasing  the  temperature 
of  the  top  portion  of  the  pavement  which  decreases  its  shear  strength  and  thus 
facilitates  easier  removal.  The  material  is  then  either  taken  to  a  batch  plant 
or  recycled  in-place.  Additional  asphalt,  rejuvenators ,  or  virgin  materials  are 
added  before  the  reworked  asphaltic  concrete  is  laid  down  and  compacted. 

During  the  past  several  decades  asphaltic  concrete  pavements  have  been  re¬ 
cycled  or  improved  using  heater  planer  and  heater  scarifier  techniques  with 
little  knowledge  of  the  temporal  and  spatial  temperature  profile  created  during 
the  process.  This  information  is  essential  if  such  methods  are  to  be  optimized 
with  respect  to  constraints  on  the  maximum  allowable  temperature  induced  in  the 
asphaltic  concrete  and  the  energy  (usually  fossil  fuel)  consumed  during  the 
process. 

Few  published  studies  have  been  conducted  on  the  temporal-spatial  tempera¬ 
ture  profiles  resulting  from  the  pavement  heating  process,  and  the  literature 
contains  little  documented  information  on  the  thermal  properties  of  asphaltic 
concrete  pavement.  The  heating  process  must  produce  pavement  temperatures  that 
fall  within  certain  constraints.  Pavement  temperature  limits  of  interest  are: 
(1)  the  temperature  at  the  surface  must  be  less  than  the  flash  point  of  asphalt 
and;  (2)  the  temperature  at  the  depth  of  removal  must  be  sufficiently  high  for 
easy  removal.  A  mathematical  thermal  model  of  asphaltic  concrete  pavement  is 
necessary  to  predict  the  temporal-spatial  temperature  profiles  that  result  from 
the  heating  process.  The  impetus  for  studying  asphaltic  pavement  thermal  pro- 


perties  stems  from  the  lack  of  studies  concerning  how  the  thermal  properties  of 
asphaltic  pavement  vary  with  different  mix  designs  and  aggregate  types.  The 
objectives  of  this  project  were  to  model  the  heating  process  and  study  the 
thermal  properties  of  asphaltic  pavements.  Four  steps  were  taken  to  achieve 
these  objectives: 

1)  Obtain  in  the  laboratory  an  expected  range  of  asphaltic  pavement  ther¬ 
mal  properties  and  determine  how  the  thermal  properties  vary  with 
asphalt  content,  gradation,  and  aggregate  type  (Part  I  of  this  report). 

2)  Develop  a  heat  transfer  model  and  computer  code  specifically  for  the 
type  of  heating  process  encountered  in  in-place  surface  recycling 
methods  (Part  II  of  this  report). 

3)  Validate  the  model  by  comparing  the  temporal-spatial  temperature  data 
observed  in  field  tests  to  those  predicted  data  determined  from  the 
computer  program  (Part  III  of  this  report). 

4)  Conduct  a  parametric  study  of  the  heating  process  to  determined  how 
sensitive  the  process  is  to  the  variation  of  thermal  and  physical 
parameters  (Part  IV  of  this  report). 

Part  V  of  this  report  is  a  summary  of  important  results  and  findings  of  the 
study. 
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PART  I 


THERMAL  PROPERTIES  OF 
SOME  ASPHALTIC  PAVEMENTS 


AND  MIXES 


CHAPTER  I-A  INTRODUCTION 


Although  thermal  properties  of  asphalt  and  asphaltic  concrete  have  been 
reported  in  the  literature  for  at  least  55  years  (14),  much  of  the  data  is  un¬ 
documented  in  that  variables  that  can  be  expected  to  influence  the  thermal  pro¬ 
perties  such  as  asphalt  content,  density,  and  aggregate  type  and  gradation  were 
not  published.  Prediction  of  in  situ  temporal-spatial  temperature  profiles 
created  by  surface  heating  requires  that  thermal  properties  of  the  asphaltic 
concrete  be  known.  Furthermore,  the  expected  range  of  these  properties  for  a 
wide  variety  of  asphaltic  concrete  mixes  would  also  be  valuable. 

This  part  of  the  report  reviews  the  literature  relevant  to  pavement  thermal 
properties,  describes  an  apparatus  and  procedure  by  which  thermal  properties 
were  measured,  and  reports  the  results  of  experiments  carried  out  on  six  aggre¬ 
gate  mixes  prepared  in  the  laboratory  and  on  field  samples  taken  from  Myrtle 
Beach  Air  Force  Base  and  Loring  Air  Force  Base. 

CHAPTER  I-B  REVIEW  OF  THE  LITERATURE 

A  review  of  the  literature  concerning  asphaltic  pavement  thermal  properties 
revealed  few  studies  conducted  solely  on  the  thermal  properties  of  asphaltic 
pavements.  However,  three  allied  areas  were  found  that  use  the  asphaltic  pave¬ 
ments  thermal  properties  as  parameters  in  temperature-predicting  thermal  models. 
The  allied  areas,  which  are  reviewed  below,  are  (1)  effects  of  climate  on  pave¬ 
ment  temperatures,  (2)  asphaltic  pavement  cooling  rates,  and  (3)  thermal  models 
of  in-place  surface  recycling. 

Effects  of  Climate  on  Pavement  Temperatures 

The  effects  of  climate  on  pavement  temperatures  can  be  subdivided  into  frost 
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penetration  studies  and  pavement  temperature  studies. 


Frost  Penetration  Research.  In  cold  climates,  the  deteriorating  effects  of 
frost  action  under  pavements  become  increasingly  important  as  maintenance  costs, 
traffic  frequency,  and  wheel  loadings  increase.  Research  has  been  conducted  to 
offer  a  better  understanding  of  the  processes  involved  in  frost  penetration. 

Carlson  and  Kersten  (7)  applied  heat  transfer  theory  to  the  pavement  and 
soil  subgrade  to  develop  a  model  used  to  predict  frost  penetration  below  as¬ 
phaltic  pavement.  Input  parameters  include  the  thermal  conductivity  of  the 
pavement,  the  latent  heat  of  fusion  of  water,  and  the  surface  freezing  index  of 
the  pavement.  The  surface  freezing  index  is  defined  as  the  time-average 
difference  in  temperature  between  the  ground  surface  and  the  freezing  point 
during  the  freezing  period.  The  value  of  the  thermal  conductivity  (a  measure 
of  the  ease  which  heat  energy  will  flow  through  a  medium)  used  by  Carlson  and 
Kersten  was  1.44  w/m°C  (obtained  from  thermal  property  tests  made  at  the  Univer¬ 
sity  of  Minnesota  in  1948). 

Aldrich  (2)  followed  up  on  Carlson  and  Kersten' s  work  and  developed  an  im¬ 
proved  method  of  calculating  frost  penetration.  His  method  includes  the  effects 
of  specific  heat  (the  amount  of  energy  necessary  to  change  the  temperature  of  a 
unit  volume  of  a  substance  by  one  degree).  Aldrich  used  a  thermal  conductivity 
value  of  1.45  w/m°C  in  his  analysis. 

Dempsey  and  Thompson  (15)  used  a  one-dimensional  transient  heat  flow  model 
to  predict  temperature  with  depth  in  conjunction  with  highway  frost  studies. 
They  reported  that  the  accuracy  of  the  temperatures  predicted  by  their  model 
depended  more  on  the  quality  of  the  input  data  than  the  numerical  method  of 
calculation.  They  also  pointed  out  the  importance  of  accurately  defining  the 
boundary  condition  at  the  pavement  surface,  because  it  is  this  input  which  is 
the  major  factor  contributing  to  the  heat  transfer  process. 
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Minimum  and  Maximum  Pavement  Temperatures  Due  to  Weather.  Because  the 
strength  and  stability  of  an  asphaltic  pavement  are  related  to  its  temperature, 
research  has  been  conducted  to  determine  the  minimum  and  maximum  temperatures  of 
pavement  in  natural  field  conditions.  Studies  indicate  that  given  the  thermal 
properties  of  the  asphaltic  concrete  pavement  and  the  history  of  the  weather  in 
the  area,  pavement  temperatures  can  be  predicted  by  heat  transfer  models. 

Barber  (3)  used  the  diffusion  equation  to  predict  in-situ  pavement  tempera¬ 
tures.  Weather  parameters  involved  in  the  diffusion  equation  include  wind  speed, 
precipitation,  air  temperature,  and  solar  radiation.  Barber  averaged  several 
published  values  of  asphaltic  pavement  thermal  conductivity  to  obtain  1.21  w/m°C. 

Southgate  and  Been  (25)  analyzed  the  pavement  temperatures  recorded  at 
various  depths  in  a  0.3  m  thick  asphaltic  pavement.  They  found  that  a  fourth 
order  polynomial  fit  the  temperature  as  a  function  of  depth  data.  No  attempts 
were  made  to  analyze  the  data  using  heat  transfer  theory.  Rumney  and  Jimenez 
(23)  recorded  pavement  temperatures  with  depth  for  one  year  in  Arizona.  Using 
the  temperature  data,  they  empirically  derived  curves  which  could  be  used  to 
estimate  the  pavement  temperature  at  any  depth  knowing  the  maximum  air  tempera¬ 
ture  and  average  daily  radiation  rates. 

Straub,  Schenk,  and  Przybycien  (28)  used  weather  data  and  pavement  tempera¬ 
ture  data  gathered  from  a  test  site  for  a  period  of  one  year  to  develop  a  heat 
transfer  model  which  could  predict  pavement  temperatures.  Their  forward  differ¬ 
ence,  one-dimensional  transient  heat  flow  model  used  air  temperature  and  solar 
radiation  data  and  the  thermal  and  physical  properties  of  a  specific  locale  to 
predict  maximum  and  minimum  pavement  temperatures.  They  found  that  the  pave¬ 
ments  had  a  substantial  temperature  gradient  and  no  one  temperature  was  repre¬ 
sentative  of  that  of  the  pavement.  They  used  Barber's  averaged  value  for  thermal 
conductivity  of  1.21  w/m°C. 


Christison  and  Henderson  (11)  also  used  a  finite  difference  approximation 
to  predict  temperatures  in  asphaltic  concrete  pavement.  They  assumed  a  value  of 
1.45  w/m°C  for  the  thermal  conductivity  of  asphaltic  concrete  and  stated  that 
"the  thermal  conductivity  (k)  and  specific  heat  (c)  of  asphaltic  concrete  paving 
mixtures  vary  within  narrow  limits  and  for  practical  purposes  can  be  considered 
independent."  However,  they  did  not  indicate  the  basis  for  this  statement  nor 
their  source  for  the  assumed  thermal  properties. 


Cooling  Rates  of  Hot-Laid  Asphaltic  Pavements 

During  paving  operations,  a  hot-laid  asphaltic  pavement  must  be  compacted 
while  the  pavement  is  cool  enough  to  have  the  stability  needed  to  support  a 
roller  or  other  compacting  device  and  warm  enough  to  achieve  a  specified  density 
with  a  minimum  amount  of  compactive  effort.  Compaction  attempts  made  when  the 
pavement  is  too  cool  result  in  either  more  compactive  effort  needed  to  reach 
a  specified  density  or  the  inability  to  reach  the  specified  density  at  all. 
Strength  and  durability  problems  have  been  associated  with  improper  compaction 
of  asphaltic  pavements. 

The  amount  of  time  available  for  compaction  depends  upon  the  cooling  rate 
of  the  hot  asphalt  mat.  The  cooling  rate  depends  upon  many  thermal  and  physical 
properties  of  the  asphaltic  pavement,  the  temperature  of  the  hot  pavement  mix, 
the  temperature  of  the  surrounding  air,  the  temperature  of  the  existing  pavement, 
the  wind  speed,  and  the  amount  of  solar  radiation  striking  the  pavement  surface. 

Corlew  and  Dickson  (12)  considered  all  of  these  factors  in  their  one-dimen¬ 
sional  transient  heat  flow  equation  to  predict  the  temporal-spatial  temperature 
field  within  a  cooling  layer  of  freshly  placed  hot-mix  asphaltic  concrete.  The 
pertinent  weather  data  and  the  temperature  of  the  hot-laid  asphaltic  pavement 
were  measured  at  the  site,  while  the  thermal  and  physical  properties  were  as¬ 
sumed  using  values  reported  in  the  literature  as  a  guide.  Using  a  finite  dif- 
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ference  technique,  they  were  able  to  predict  temperatures  in  a  six  centimeter 

layer  of  asphaltic  concrete  within  seven  degrees  centigrade  of  measured  tempera-  ; 

tures  in  the  range  at  which  asphaltic  concrete  is  compacted  (greater  than  90°C). 

With  the  impetus  of  extending  the  paving  season  into  the  cooler  weather 
seasons,  Frenzel  et  al .  (16)  developed  a  computer  analysis  to  study  the  effect 
of  preheating  an  existing  asphaltic  layer  on  the  cooling  of  an  overlay  put  down 
over  the  preheated  base.  They  found,  analytical ly ,  that  preheating  the  base 
increased  the  cooling  down  period  in  thin  overlays  sufficiently  to  allow  a  com- 

i 

.  I 

paction  window  long  enough  to  make  paving  in  early  spring  and  late  fall  feasible 
Later  Corlew  and  Dickson  (13)  combined  the  base  preheat  model  developed  by 
Frenzel  et  al.  and  the  previously  described  cooling  model  developed  by  Corlew 
and  Dickson  (12)  to  predict  the  temporal-spatial  temperature  field  of  an  as¬ 
phaltic  concrete  layer  over  a  preheated  base  in  a  bench  scale  laboratory  test. 

A  sample  with  a  10  cm  diameter  base  was  heated  by  a  direct-fired  propane  heater. 

Agreement  between  predicted  and  measured  temperatures  was  considered  to  be  good. 

Assumed  thermal  properties  of  asphaltic  concrete  were  used  in  the  analysis. 

Another  study  of  asphaltic  pavement  cooling  rates  was  conducted  by  Wolfe 
and  Colony  (30)  at  the  University  of  Toledo.  They  developed  a  computer  simu¬ 
lation  method  to  predict  cooling  rates  using  the  same  weather  data  and  material 
property  variables  as  Corlew  and  Dickson.  The  thermal  properties  of  asphaltic 
pavement  were  taken  from  published  values.  Later  Wolfe,  Heath,  and  Colony  (31) 
continued  the  research  and  further  verified  the  computer  model.  In  this  study, 
the  thermal  conductivity  of  seven  different  mix  designs  of  ashphaltic  concrete 
was  measured  using  a  guarded  hot-plate  thermal  conductivity  instrument  as 
specified  in  ASTM  C- 177.  Thermal  conductivity  values  ranged  from  0.58  to  1.01 
w/m°C.  The  specific  heat  of  the  seven  mixes  was  also  measured  and  found  to 
range  from  878  to  962  J/kg°C. 
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Thermal  Models  of  In-Place  Surface  Recyclinc 


In-place  surface  recycling  of  asphaltic  pavements  involves  reworking  the 
surface  to  a  depth  of  approximately  2.5  centimeters  using  a  heater  scarifier. 
This  operation  may  involve  the  addition  of  new  materials  (or  recycled  materials) 
including  aggregates,  modifiers,  and/or  asphalt  cement.  The  reworked  material 
is  then  compacted,  and  sometimes  a  seal  coat  is  applied.  Typically,  heater 
planers  and  heater  scarifiers  developed  by  contractors  use  propane  as  a  fuel  to 
fire  a  grid  of  torches.  The  grid  is  on  the  order  of  3  m  wide  and  5  m  or  more 
long  and  is  attached  to  a  self-propelled  machine.  The  torch  grid  is  covered  on 
the  top  and  sides  to  diminish  heat  loss.  Typically,  the  machine  advances  at 
the  rate  of  about  4.5  m/min  giving  a  heat  exposure  time  on  the  order  of  one 
minute.  To  prevent  combustion  of  the  pavement,  the  temperature  at  the  pavement 
surface  should  be  limited  to  about  230°C.  Contractors  claim  that  under  these 


conditions,  temperatures  at  a  depth  of  about  2.5  cm  in  the  asphaltic  concrete 
are  sufficiently  high  to  allow  easy  scarification  and/or  removal  of  pavement  to 
this  depth.  To  increase  the  depth  of  influence  of  the  heat  source,  a  technique 
termed  "soaking"  is  sometimes  used.  This  method  involves  heating  the  surface  as 
described  above  and  then  insulating  the  pavement  surface.  This,  in  principle, 
allows  the  heat  to  "soak"  in  thus  producing  greater  heat  penetration  without 
exceeding  the  230°C  surface  temperature  requirement. 

One  of  the  few  published  studies  undertaken  in  the  area  of  heat  transfer  in 
asphaltic  concrete  recycling  seems  to  refute  many  of  the  contractor  claims  men¬ 
tioned  above.  Carmichael,  Boyer,  and  Hokanson  (8)  modeled  an  asphaltic  concrete 
pavement  as  a  semi-infinite  solid  and  used  a  forward  difference  numerical  method 
to  solve  the  governing  differential  equation.  They  assumed  that  the  thermal 
conductivity,  density,  and  specific  heat  of  asphaltic  concrete  were  independent 
of  temperature  and  did  not  vary  from  point  to  point.  The  authors,  using  real¬ 
istic  source  and  initial  pavement  temperatures  along  with  realistic  pavement 
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parameters,  found  that  even  with  an  induced  surface  temperature  of  540°C  and  a 
30  second  exposure  time,  there  was  no  increase  in  temperature  of  the  pavement  at 
a  depth  of  1.6cm  (See  Figure  1-1). 

It  was  also  shown  that  for  surface  temperatures  between  200  and  300°C,  the 
temperature  of  the  pavement  hardly  changed  from  its  initial  condition.  This  is 
illustrated  in  Figure  1-2.  Although  it  can  be  argued  that  this  study  used  re¬ 
strictive  assumptions,  particularly  with  respect  to  the  parameters  used,  and 
that  the  model  has  not  been  verified,  it  does  indicate  that  the  current  develop¬ 
ment  of  pavement  heating  techniques  and  the  present  use  of  heater  planer  and 
heater  scarifier  equipment  ignores  some  of  the  basic  constraints  of  asphaltic 
concrete  recycl i ng--namely  that  to  obtain  sufficient  temperature  increases  at 
depth  it  appears  that  the  surface  must  be  heated  to  a  temperature  greatly  in  ex¬ 
cess  of  that  at  which  damage  to  the  pavement  can  occur.  Carmichael  et  al.  point¬ 
ed  out  that  it  is  important  that  the  full  depth  of  pavement  be  properly  heated 
before  scarifying  because  "cold  asphalts  do  not  bond  well  and  any  attempt  to 
force  cold,  asphalt  coated  aggregates  together  after  scarifying  them  by  tensile 
failure  can  only  lead  to  future  problems." 

Thermal  Properties  of  Asphaltic  Concrete  and  Its  Components 

Much  of  the  literature  cited  above  used  assumed  values  of  the  thermal 
properties  of  asphaltic  concrete.  However,  as  shown  in  Tables  1-1  and  1-2 

measured  values  of  conductivity  and  specific  heat  of  asphalt  and  asphaltic  con¬ 
crete  as  well  as  minerals  and  rocks  commonly  used  as  asphaltic  concrete  aggre¬ 
gate  are  available.  Reported  values  of  thermal  conductivity  range  from  0.74  to 
2.88  w/m°C  for  asphaltic  concrete  and  from  0.14  to  0.17  w/m°C  for  "pure  asphalt." 

.7  _  7 

Diffusivity  values  ranging  from  5.8*10  to  14.4*10  m2/s  have  been  reported. 
The  range  of  reported  values  of  specific  heat,  c,  is  879  to  963  J/kg°C  for  as¬ 
phaltic  concrete. 
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Figure  1-2.  Pavement  Temperature  vs.  Depth  as  a  Function 
of  Source  Temperature  and  Exposure  Time  (8). 


Table  1-1. 

Some  Published 

Thermal  Properties 

of  Asphaltic  Concrete. 

Property 

k(w/m°C) 

a(m2/s) 

pc( J/m3°C) 

Ref. 

Comment 

1.454 

2.88 

14.4*10~7 

1. 41*106 

2. 00*106 

32 

19 

18°C ,  Dry 

2.28 

11.5*10  ' 

1. 97*106 

19 

38°C ,  Dry 

1.21 

5.75*10  ' 

8 

Obtained  from  the 
Asphalt  Institute 

0.74-0.76 

6 

20°C-56°C 

0.167-0.172 

14 

"Pure  Asphalt" 
2C°C-80°C 

0.65-0.75 

-7 

14 

Asphalt  used  in 
street  paving 

1.37-1.75 

7-8*10  ' 

c=879-963 

31 

0.80 

io*io'7 

c=879 

14 

26 

100°C 

1.21 

c=920 

3 

1.2 

5.9*10 

2. 07*106  12, 

13 

80°C-149°C 

c=921 

1.5 

2 

0.82-2.32 

22 

0.14-0.17 

0=1528-2561 

24 

Asphaltic  Bitumen 
free  of  paraffin 
0°C-300°C 
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Table  1-2.  Some  Published  Thermal  Properties  of  Asphaltic 
Concrete  Aggregate. 


Property 


k(w/m°C) 

a(m2/s) 

c(J/kg°C) 

Ref. 

Comment 

790;  1000 

18 

Cal  cite  0°C;  200°C 

930 

18 

Dolomite  60°C 

698;  969 

18 

Quartz  0°C;  200°C 

1000 

18 

Limestone  58°C 

680 

18 

Limestone,  mean  of 

3  at  50°C 

830 

18 

Limestone,  mean  of 

10  at  65°C 

700;  970 

18 

Quartzite  0°C; 

200°C 

650;  950 

18 

Granite  0°C;  200°C 

850;  1040 

18 

Basalt  0°C ;  200°C 

930 

18 

Sandstone  59°C 

700;  870 

18 

Diabase  0°C;  200°C 

710;  1000 

18 

Slate  0°C ;  200°C 

1.82 

17 

Average  value  for 

sand  and  gravel 

2.86 

5 

Cal  cite  100°C 

6.45 

5 

Quartz  100°C 

2.37 

5 

Granite  100°C 

1.8-2. 2 

5 

Basalt 

2. 0-3. 4 

5 

"Compact" 

Limestone 

1.1-2. 2 

5 

"Porous" 

Limestone 

1.8 

5 

Slate  100°C 

3.99 

5 

Dolomite  100°C 

5.2 

5 

Quartzite  100°C 

1.8-2. 8 

5 

Granite-Gneiss 

2.7 

5 

Grani te-Schi st 

2. 6-4. 5 

5 

Hard  Sandstone 

2.1 

1.1 

2*10"7 

5 

5 

Diabase  100°C 
Quartz  Sand 

9*10 

5 

Sandy  Soil 
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The  predominant  component  of  asphaltic  concrete  is  aggregate.  The  thermal 
properties  of  the  particular  aggregate  used  in  a  given  asphaltic  concrete  can  be 
expected  to  have  an  important  effect  on  the  properties  of  the  mix.  The  liter¬ 
ature  indicates  that  common  minerals  and  rocks  used  as  aggregate  can  have  con¬ 
ductivities  which  range  from  about  1  to  8  w/m°C  (Table  1-2);  the  range  of  c  can 
be  expected  to  be  on  the  order  of  650  to  1040  J/kg°C.  Wolfe  et  al.  (31)  suggest 
920  J/kg°C  as  an  average  value  based  on  their  tests.  Limited  data  on  the  diffu- 
sivity  of  the  fine  aggregate  component  indicate  a  is  on  the  order  of  2  to  9*10  7 
m2/s. 

Thermal  properties  of  asphaltic  concrete  reported  in  the  literature  do  not 
indicate  the  mineralogy,  grain  size  distribution,  or  density  of  the  aggregate 
used  in  the  mix;  nor  have  the  fraction  of  asphalt  of  density  of  mix  been  report¬ 
ed.  Certainly  these  factors  can  be  expected  to  influence  the  resulting  thermal 
porperties  of  asphaltic  concrete  and  thus  the  transfer  of  heat  through  and  as¬ 
phaltic  concrete  pavement. 

CHAPTER  I-C  MEASUREMENT  OF  THERMAL  PROPERTIES 
Properties  Required 

The  transient  heat  conduction  process  in  a  solid  body  is  described  by  the 
energy  balance  equation, 

Rate  of  Rate  of  Heat  Rate  of  Rate  of  Change 

Heat  In-  +  Generation  by  =  Heat  Out-  +  In  Internal  (1) 

Flow  Internal  Sources  Flow  Energy 


Without  any  internal  heat  generation  in  the  solid,  this  equation  can  be  written 
algebraically  in  rectangular  coordinates  as 


where  x,  y,  z  are  spatial  coordinates  (m) 
T  =  temperature  (°C) 
k  =  thermal  conductivity  (w/m°C) 
p  =  density  (kg/m3) 
c  =  specific  heat  (J/kg°C) 
t  =  time  (s) 


Field  heaters  are  typically  about  5  m  long  by  3  m  wide.  Since  these  dimensions 
are  much  greater  than  the  thickness  of  the  pavement,  the  conduction  heat  trans¬ 
fer  can  be  assumed  to  be  one-dimensional.  Therefore,  with  the  z-direction  being 
vertical , 


3T 

3x 


(3) 


and  Equation  (2)  becomes, 
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by  defining  the  thermal  diffusivity,  a,  as  a  measure  of  the  rate  at  which  a  tem¬ 
perature  change  spreads  through  a  solid  or 

a  =  k  /  p  c  (5) 

and  kz  =  k.  Equation  (4)  reduces  to 
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This  equation  is  the  one-dimensional  transient  heat  conduction  equation, 
and  in  similar  forms  can  be  used  in  other  transient  energy  flow  situations. 
One  example  of  transient  energy  flow  familiar  to  geotechnical  engineers  is 
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Terzaghi's  one-dimensional  consolidation  theory.  In  Equation  (6)  the  thermal 


diffusivity,  a,  is  analogous  to  the  coefficient  of  consolidation,  c^,  and  is  an 
indicator  of  the  rate  of  temperature  (pore  water  pressure)  dissipation  within  a 
solid  (soil  stratum). 

Equation  (6)  can  be  written  in  finite  difference  form  so  that  the  temporal - 
spatial  temperature  field  can  be  calculated  for  a  pavement  system  having  several 
layers,  provided  the  initial  and  boundary  conditions  and  the  thermal  properties, 
p,  c,  and  k  (or  a)  are  known  for  each  layer.  The  mathematics  and  computer  cod¬ 
ing  for  predicting  the  temperature  profile  with  time  for  a  layered  asphaltic  con¬ 
crete  pavement  system  for  a  known  heat  input  are  routine  provided  the  asphaltic 
concrete  thermal  properties  are  known.  This  chapter  describes  the  laboratory 
apparatus  and  procedure  used  to  measure  three  thermal  properties:  thermal  dif¬ 
fusivity,  thermal  conductivity,  and  specific  heat. 

Apparatus  and  Procedure 

Many  methods  of  measuring  thermal  properties  have  been  used  over  the  years. 
Most  methods  are  for  measuring  a  single  property  at  a  time,  and  they  usually 
utilize  closed  form  solutions  for  steady-state  and  transient  cases.  The  method 
chosen  for  this  research  has  some  special  advantages. 

Beck  and  Al-Araji  (4)  developed  a  method  of  thermal  property  measurement 
capable  of  measuring  the  thermal  conductivity,  thermal  diffusivity,  and  specific 
heat  in  a  single  test  cycle.  The  thermal  conductivity  is  considered  constant  in 
this  method.  The  method  requires  the  integration  of  some  temperature  data  which 
can  be  easily  accomplished  with  a  microcomputer.  The  test  procedure  is  quick, 
simple,  and  requires  minimal  equipment. 

In  concept,  this  method  involves  applying  a  given  amount  of  heat,  Q,  to  the 
top  surface  of  a  flat  specimen  of  known  thickness,  and  then  measuring  tempera¬ 
tures  at  the  top  (T^)  and  bottom  (T^)  surfaces  of  the  insulated  specimen  over 
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time  (see  Figure  1-3).  The  heat  is  applied  with  a  calorimeter--an  insulated 
copper  disk  of  well  known  physical  and  thermal  properties  and  at  a  uniform  tem¬ 
perature  different  than  that  of  the  specimen  (Figure  I-3a).  When  the  calori¬ 
meter  and  specimen  are  brought  into  contact,  the  top  of  the  specimen  will 
immediately  assume  the  temperature  of  the  calorimeter,  and  then  will  begin  to 
cool  as  the  heat  is  transferred  to  the  specimen  (Figure  I  -  3C . ) .  At  the  same 
time,  the  temperature  at  the  bottom  of  the  specimen  will  begin  to  rise.  Pro¬ 
vided  there  are  no  heat  losses,  the  thermal  energy  will  continue  to  travel  into 
the  specimen  until  the  specimen  and  the  calorimeter  reach  the  same  temperature, 
called  which  occurs  at  time  =  t^.  Figure  1-4  is  an  illustration  of  how  the 
temperatures  T  and  are  a  function  of  time  and  how  at  time  =  t^,  =  Tb  =  T^. 

Beck  and  Al-Araji  (4)  showed  that  for  flat-plate  geometry,  the  thermal 
properties  could  be  calculated  from  the  following  equations: 


a 


L2<Tf-  V 

t=tf 

;  f(T  -  t  )  dt 
t=o  *  D 


(7) 


k  = 


t=t. 

I 

t=0 


IS. 


<Tt- 


V 


dt 


c  =  - 3 - 

P  L  (Tf  -  T. ) 

where  a  =  thermal  diffusivity 

L  =  thickness  of  asphaltic  pavement  disk 
=  final  pavement  temperature 


(8) 


(9) 


T..  =  initial  pavement  temperature 


Figure  3a.  Calorimeter. 


insulation. 

Copper  disk  with  known  density 
k,  c.  and  thickness  (L  ). 


Figure  3b.  Sample  holder. 


Top  and  bottom  thermocouples 
measured  every  time,  dt. 


Pavement  disk  with  known 
density  and  thickness  (L). 

Insulation. 


Figure  3c. 


Calorimeter  in  place 
atop  the  sample  holder. 


Figure  1-3 (a,  b,  and  c).  Schematic  of  Calorimeter ,  S ample  Holder 

and  How  the  Two  Fit  Together. 


time 
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T^,  =  temperature  at  the  top  and  bottom  of  the  pavement 


sample,  respectively,  at  time  t 
dt  =  specified  time  increment  (15  seconds) 

Q  =  amount  of  heat  transferred  from  the  calorimeter 
to  the  asphaltic  pavement  sample 
c  =  specific  heat 

p  =  previously  (independently)  measured  density  of  the 
asphaltic  pavement  sample  (ASTM  D1188) 
k  =  thermal  conductivity  of  the  asphaltic  pavement  sample 

To  find  Q,  the  amount  of  heat  transferred  from  the  calorimeter  to  the  sam¬ 
ple,  note  that  all  surfaces  are  insulated  when  the  calorimeter  and  sample  are 
brought  into  contact  (see  Figure  I-3c.).  Therefore,  all  the  thermal  energy 
available  in  the  calorimeter  will  be  transferred  to  the  pavement  sample  until 
the  calorimeter  and  the  sample  reach  the  same  temperature  (see  Figure  1-4).  The 
amount  of  heat  transferred,  Q,  can  be  calculated  using  the  following  equation: 

0  -  pcccLc(Wc  (10) 

where  the  subscript  refers  to  the  calorimeter  and  cc  refers  to  the  specific 
heat  of  the  calorimeter  at  the  average  temperature.  The  procedure  can  be  car¬ 
ried  out  in  five  steps: 

1.  Bring  into  contact  an  insulated  calorimeter  of  known  thickness,  tempera¬ 
ture,  and  thermal  properties  with  an  insulated  specimen  of  known  thick¬ 
ness,  temperature,  and  density. 

2.  Measure  the  temperature  at  the  top  and  the  bottom  of  the  specimen  from 
the  time  of  contact  (time  =  0)  until  all  the  thermal  energy  available 
in  the  calorimeter  is  transferred  into  the  specimen  (time  =  t^.  when  T^ 
=  T^  =  T^,  see  Figure  1-4). 
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3.  Evaluate  the  integral  in  Equations  6  and  7  by  summing  the  temperature 
difference,  -  T^,  from  time  =  0  time  =  t^  and  multiplying  this  sum 
by  the  sampling  interval,  dt,  to  get: 


4.  Calculate  Q,  the  amount  of  heat  energy  transferred  from  the  calorimeter 
to  the  specimen  (see  Equation  10). 

5.  Calculate  the  thermal  diffusivity  (a),  thermal  conductivity  (k),  and 
specific  heat  (c)  using  Equations  7,  8,  and  9,  respectively. 

The  test  apparatus  used  in  this  study  consists  of  an  insulated  copper 
calorimeter,  a  calorimeter  heater  with  heater  controller,  an  insulated  pavement 
sample  holder,  and  a  temperature  measuring  system  with  five  copper-constantan 
thermocouples.  The  calorimeter  (shown  in  Figure  1-5)  is  composed  of  a  pure 
copper  disk  (having  accurately  known  thermal  properties)  10.16  cm  in  diameter  by 
3.18  cm  thick  surrounded  by  1.2  cm  of  firebrick  insulation  on  the  perimeter  and 
on  one  face.  An  oak  wood  holder  surrounds  the  insulation.  A  single  thermocou¬ 
ple  placed  in  the  center  of  the  copper  disk  is  used  to  monitor  and  control  the 
temperature  of  the  calorimeter  with  the  heater.  During  the  calorimeter  heating 
process,  the  heater  controller  regulates  the  power  to  the  calorimeter  heater 
(Figure  1-6)  in  such  a  way  as  to  bring  the  calorimeter  to  a  constant  pre¬ 
scribed  temperature. 

The  10.16  cm  diameter  asphaltic  concrete  pavement  sample  is  held  in  a 
holder  which  insulates  the  perimeter  and  the  bottom  surface  of  the  sample  (Fig¬ 
ure  1-7).  Two  copper-constantan  thermocouples  are  placed  on  both  the  top  and 
bottom  surfaces  of  the  pavement  sample  with  the  top  thermocouples  recessed  in 
shallow  grooves  sawn  in  the  top  of  the  pavement  sample.  These  thermocouples 
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Figure  1-5.  Calorimeter. 


Figure  1-6.  Calorimeter  Heater. 
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lead  to  a  temperature  measurement  system  made  by  Interactive  Microware  (see 
Appendix  I-A  and  Ref.  29)  where  the  two  top  thermocouples  are  averaged  giving 
Tj.,  and  the  two  bottom  thermocouples  are  averaged  giving  at  each  specified 
time  interval,  dt. 

The  TEMPSENSE  Temperature  Monitoring  System  (29)  in  conjunction  with  an 
Apple  11+  computer  has  user-controlled  software  which  facilitates  the  calibra¬ 
tion  of  the  thermocouples  and  allows  the  operator  to  select  the  sampling  fre¬ 
quency  and  the  total  duration  of  the  temperature  measurements.  The  time  of  each 
measurement  is  displayed  on  the  screen  along  with  the  temperature  of  each  thermo¬ 
couple  measured.  The  data  are  automatically  stored  on  a  floppy  disk  for  later 
analysis. 

To  run  a  test,  the  calorimeter  is  placed  on  the  heater,  and  the  desired 
temperature  of  the  calorimeter  is  set  on  the  heater  controller.  The  thickness, 
L,  of  the  previously  prepared  pavement  sample  disk  is  measured,  and  the  bottom 
thermocouples  are  fastened  onto  the  center  of  the  bottom  surface  of  the  sample 
with  a  viscid,  high-conductivity  silicon  paste  (Omegatherm  201  from  Omega  Engi¬ 
neering,  Inc.)  and  masking  tape.  The  sample  is  then  placed  in  the  insulated 
holder  and  a  thin  film  of  the  silicon  paste  is  applied  to  the  top  surface.  The 
top  thermocouples  are  then  laid  into  the  recessing  grooves,  and  an  insulating 
collar  is  placed  over  the  sample  and  thermocouples.  When  the  calorimeter  reach¬ 
es  the  prescribed  temperature,  the  pavement  sample  is  ready  for  contact  with  the 
heated  calorimeter.  The  prescribed  calorimeter  temperature  (145°C)  was  chosen 
to  be  approximately  halfway  between  the  minimum  desirable  pavement  temperature 
for  compaction  (90°C)  and  the  typical  flash  point  temperature  of  asphalt  cement 
(about  230°C) . 

The  TEMPSENSE  program  prompts  for  the  desired  sampling  interval,  the  dura¬ 
tion  of  the  temperature  measurements,  and  the  file  name  of  the  data  to  be  gath¬ 
ered.  A  fifteen  second  sampling  interval  was  chosen  to  provide  reasonable 
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precision,  and  it  was  found  that  test  durations  of  around  30  minutes  (t^  = 
min.)  were  sufficient  for  =  T^.  The  temperature  measurements  are  started,  and 
two  or  three  measurements  are  recorded  to  obtain  the  initial  temperature  of  the 
sample  (T. ).  Then  the  calorimeter  is  removed  from  the  heater  and  placed  on  the 
top  of  the  sample.  Measurements  are  made  every  15  seconds  until  the  tempera¬ 
ture  at  the  bottom  of  the  pavement  sample  becomes  equal  or  very  close  to  the 
temperature  at  the  top  of  the  pavement  sample,  i.e.  until  tf  has  elapsed.  The 
temperature  data  is  then  stored  on  a  floppy  disk,  and  the  calorimeter  is  removed 
and  placed  back  on  the  heater  in  preparation  for  the  next  test. 

The  integral  expression  in  the  denominator  of  Equations  7  and  8  is  approxi¬ 
mated  by  summing  the  temperature  difference  between  the  top  and  bottom  of  the 
sample  for  each  dt  increment  of  time  (dt  =  15  seconds)  over  the  duration  of  the 
test  (see  Equation  11).  The  duration  of  the  test  is  defined  as  the  time  that 
elapses  from  the  time  the  calorimeter  and  sample  first  contact  (time  =  0)  to  the 
time  an  equilibrium  temperature  is  reached  throughout  the  calorimeter  and  sample 
(time  =  tf  and  Tt  =  Tfa). 

By  using  the  collected  data  and  the  previously  measured  density  of  the  pave 
ment  sample,  Equations  7,  8,  and  9  are  solved  for  the  thermal  diffusivity,  ther¬ 
mal  conductivity,  and  the  specific  heat,  respectively. 

CHAPTER  I-D  PREPARATION  OF  LABORATORY  ASPHALTIC  CONCRETE  SAMPLES 

This  chapter  provides  a  brief  description  of  asphaltic  concrete  and  the 
variables  considered  in  asphaltic  concrete  mix  design.  The  reasons  for  choosing 
the  specific  mix  designs  used  in  this  study  are  explained,  and  the  designs  are 
then  summarized  in  Table  1-3  and  Figure  1-8. 

Typically,  hot-laid  asphaltic  concrete  pavements  contain  approximately  5% 
by  weight  asphalt  cement  and  95%  by  weight  graded  aggregate.  Each  aggregate 
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Table  1-3.  Summary  of  Mix  Designs  (each  prepared  with  3.5%,  5.0%, 
and  6.5%  asphalt  content). 


Aggregate  Type 

Gradations 

Dense 

Open 

Base 

Limestone 

V 

V 

V 

Granite-Gneiss 

V 

Expanded  Shale 

V 
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Figure  1-8.  Grain  Size  Distributions  for 


type  and  gradation  has  associated  with  it  a  percentage  of  asphalt  cement  that 
provides  an  optimum  mixture  of  durability,  stability,  and  workablity. 

The  range  of  the  variables,  aggregate  type,  aggregate  gradation,  and  per¬ 
cent  asphalt  cement  used  for  this  study  were  chosen  for  two  purposes,  (1)  to 
cover  a  wide  range  of  "typical"  pavement  mixes,  and  (2)  to  study  the  effects  of 
asphalt  content,  particle  size/shape  distribution,  aggregate  mineralogy,  and  den¬ 
sity  on  the  thermal  properties  of  asphaltic  concrete  pavement. 

The  type  of  aggregate  used  in  an  asphaltic  concrete  pavement,  in  a  given 
geographical  area,  depends  largely  upon  what  is  available  in  that  region.  These 
different  types  of  aggregates,  when  classified  by  their  specific  gravity,  are 
commonly  catergorized  as  lightweight,  medium  weight,  and  heavy  weight  aggregates. 
Three  types  of  aggregate  were  used  in  this  research,  (1)  a  lightweight, 
expanded  shale  aggregate,  (2)  a  medium  weight  limestone  aggregate,  and  (3)  a 
medium  weight  granite-gneiss  aggregate. 

The  aggregate  gradation  of  an  asphaltic  pavement  mix  depends  mainly  upon 
the  function  of  the  mix.  Usually,  a  very  dense  gradation  is  desirable  for  sur¬ 
face  and  base  courses  when  a  waterproof  pavement  is  needed.  Sometimes  an  open 
graded,  porous  gradation  is  used  as  the  surface  course  to  prevent  the  ponding  of 
water  on  the  surface.  Gradations  for  three  types  of  asphaltic  pavements  were 
used  in  this  study:  dense  graded  surface  course,  dense  graded  base  course,  and 
open  graded  surface  course. 

The  particle  size  distributions  (by  weight)  of  the  asphaltic  concrete  mixes 
are  shown  in  Figure  1-8.  The  base  course  and  dense  graded  surface  course  mixes 
were  similar  except  that  the  base  course  mix  has  particle  sizes  up  to  38  mm 
whereas  the  largest  particles  present  in  the  dense  graded  surface  course  mix 
were  12  mm.  The  open  graded  surface  course  mix  has  the  same  range  of  particle 
sizes  as  the  dense  graded  surface  course  mix,  but  as  Figure  1-8  illustrates, 
aggregate  size  is  much  more  uniform.  The  dense  graded  surface  course  made  with 
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the  lightweight  aggregate  shows  a  particle  size  distribution  different  from  the 
other  dense  graded  surface  courses.  This  is  because  the  curve  is  based  on  the 
percent  passing  by  weight,  and  the  lightweight  mix  contains  a  combination  of 
lightweight,  expanded  shale  aggregate  and  other  medium  weight  aggregates. 

The  percent  of  asphalt  cement  used  in  a  particular  aggregate  mix  varies 
with  the  aggregate  type  and  gradation.  No  attempts  were  made  to  optimize  the 
percentage  of  asphalt  cement  for  the  aggregate  types  and  gradations  used  in  this 
research.  Rather,  three  percentages  by  weight  of  asphalt  cement  were  used  with 
each  aggregate  mix  so  that  the  effects  of  asphalt  content  on  the  thermal  proper¬ 
ties  could  be  studied  independently  of  aggregate  type  and  gradation.  Type  AC-20 
asphalt  cement  was  used  in  all  the  mixes. 

Asphaltic  concrete  samples  were  compacted  as  specified  by  the  Marshall 
Method  of  Mix  Design  (ASTM  D1559).  Samples  of  each  mix  were  prepared  at  asphalt 
contents  of  3.5%,  5.0%,  and  6.5%.  The  bulk  specific  gravity  (density)  of  each 
sample  was  determined  according  to  ASTM  D1188.  Pavement  disks,  about  2.5  cm 
thick  were  cut  with  a  diamond  saw  from  the  10.2  cm  diameter  by  nominal  6  cm 
thick  Marshall  samples,  and  thermocouple  recessing  grooves  were  sawn  into  one 
surface  of  each  disk.  The  thermal  properties  of  each  disk  were  then  obtained 
by  the  procedure  described  previously. 

CHAPTER  I-E  PRESENTATION  AND  DISCUSSION  OF  RESULTS 

The  thermal  properties  (a,  k,  and  c)  were  obtained  for  each  asphaltic  con¬ 
crete  mix  (mixes  summarized  in  Table  1-3)  using  the  method  described  in  Chapter 
I-C.  This  chapter  presents  the  results  of  these  tests  and  indicates  the  vari¬ 
ation  of  thermal  properties  with  aggregate  type,  aggregate  gradation,  and  as¬ 
phalt  content. 

The  density-asphalt  content  relationship  for  each  of  the  mixes  is  shown  in 
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Figure  1-9.  Each  point  on  the  curves  represents  the  average  of  the  results  of 
either  four  or  five  samples.  Thermal  properties  of  each  mix  were  obtained  from 
the  disk  samples  and  the  results  for  the  thermal  conductivity  (k),  specific  heat 
(c),  and  thermal  diffusivity  (a)  are  shown  in  Figures  1-10,  I  - 11 ,  and  1-12, 
respectively,  and  in  Table  1-4.  A  comparison  of  Figures  1-8  and  1-9  indicates 
that  thermal  conductivity  for  the  six  mixes  does  not  vary  systematically  with 
density  or  asphalt  content. 

For  the  dense  graded  surface  course  mix  with  limestone  aggregate,  the  ther¬ 
mal  conductivity  increased  as  the  asphalt  content  increased  from  3.5%  to  5.0%. 
That  the  density  of  the  mix  increased  as  well  (Figure  1-9)  indicates  that  as¬ 
phalt  was  replacing  air  in  the  mix  voids.  Since  the  thermal  conductivity  of 
asphalt  is  greater  than  that  of  air,  the  conductivity  increased  slightly  (by 
about  4%).  As  the  asphalt  content  of  the  mix  was  increased  from  5.0%  to  6.5%, 
the  density  remained  constant.  These  data  suggest  that  in  this  range  of  asphalt 
content  for  this  mix,  asphalt  had  filled  the  voids  and  some  of  the  mineral 
aggregate  itself  was  being  replaced  by  asphalt.  Since  the  conductivity  of  as¬ 
phalt  is  less  than  that  of  limestone  (Tables  1-1  and  1-2),  the  thermal  conduc¬ 
tivity  of  the  mix  decreased  as  shown  in  Figure  1-10. 

The  dense  graded  surface  course  mix  with  granite-gneiss  aggregate  behaved 
very  similarly  to  the  dense  graded  surface  course  mix  with  limestone  aggregate 
as  the  asphalt  content  was  increased  from  3.5%  to  6.5%.  The  density  (Figure 
1-9)  and  the  thermal  conductivity  (Figure  1-10)  both  increased  as  asphalt  con¬ 
tent  changed  from  3.5%  to  5.0%  indicating  that  the  air  voids  in  the  mix  were 
being  filled  with  asphalt.  From  an  asphalt  content  of  5.0%  to  6.5%,  the  density 
decreased  slightly  (0.4%)  indicating  that  the  asphalt  was  displacing  the  mineral 
aggregate  thus  decreasing  the  thermal  conductivity  of  the  mix  (Figure  I - 10 ) . 
The  values  of  the  thermal  conductivity  for  the  granite-gneiss  dense  graded  sur¬ 
face  course  mixes  were  all  about  33%  lower  than  that  of  the  limestone  dense 
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Figure  1-9.  Bulk  Density  as  a  Function  of  Asphalt  Content  for 
the  Six  Aggregate  Mixes. 
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Figure  1-10.  Thermal  Conductivity  as  a  Function  of  Asphalt 
Content  for  the  Six  Aggregate  Mixes. 
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Figure  I -1 1 .  Specific  Heat  as  a  Function  of  Asphalt  Content 
for  the  Six  Aggregate  Mixes. 


DIFFUSIVITV  (m2/s )*10 


8.5*1 


3.5  5.0  6.5 


ASPHALT  CONTENT  {%  BY  WEIGHT) 


Legend 

O  BC  -  Base  Course 
A  OGSC  -  Open  Graded  Surface  Course 
□  DGSC  -  Dense  Graded  Surface  Course 

- - Limestone  Aggregate 

—  —  —Lightweight  Aggregate 
.  Granite-gneiss  Aggregate 

Figure  1-12.  Thermal  Oiffusivity  as  a  Function  of  Asphalt 
Content  for  the  Six  Aggregate  Mixes. 
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Table  1-4.  Summary  of  Physical  and  Thermal  Properties  of  the 
Six  Aggregate  Mixes. 


PROPERTIES 

(average  of  4  or  5  samples) 

Aggregate  type 
and  Gradation 

%  Asphalt  p 

(Mg/m3) 

k 

(w/m°C) 

c 

(J/kg°C) 

a 

(m2/s*10~7 ) 

Limestone 

DGSC 

3.5 

2.41 

1.53 

1040 

6. 1 

5.0 

2.47 

1.59 

968 

6.7 

6.5 

2.47 

1.47 

970 

6.2 

OGSC 

3.5 

2.10 

1.43 

1135 

6.1 

5.0 

2. 12 

1.36 

1120 

6.2 

6.5 

2.23 

1.42 

1145 

6.2 

BC 

3.5 

2.50 

1.48 

900 

6.6 

5.0 

2.51 

1.55 

802 

7.8 

6.5 

2.48 

1.66 

900 

7.3 

Granite-gneiss 

DGSC 

3.5 

2.37 

0.97 

682 

6.0 

5.0 

2.46 

1.06 

828 

5.3 

6.5 

2.45 

1.00 

677 

6.0 

OGSC 

3.  5 

2.21 

0.95 

908 

4.7 

5.0 

2.25 

0.92 

907 

4.5 

6.5 

2.28 

0.98 

900 

4.6 

Expanded  Shale 

DGSC 

3.  5 

1.58 

0.90 

1550 

3.7 

5.0 

1.62 

0.88 

1480 

3.7 

6.5 

1.63 

0.94 

1600 

3.7 

LEGEND 

DGSC  -  Dense  Graded 
OGSC  -  Open  Graded 
BC  -  Base  Course 

Surface  Course 

Surface  Course 
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graded  surface  course  mixes.  This  is  to  be  expected  because  (1)  as  shown  in 
Table  1-2,  the  thermal  conductivity  of  granite-gneiss  is  less  than  that  of  lime¬ 
stone  (10%  to  18%  less),  and  (2)  the  density  of  the  granite-gneiss  mix  is  slight¬ 
ly  (1%)  less  than  that  of  the  limestone  mix  (Figure  1-9). 

The  thermal  conductivity  of  the  limestone  base  course  mix  increased  with 
asphalt  content  over  the  entire  range  of  asphalt  content  (3.5%  to  6.5%)  at  which 
samples  were  tested.  This  behavior  was  observed  despite  the  fact  that  the  den¬ 
sity  of  the  base  course  mix  changed  little  with  asphalt  content.  It  is  believed 
that  this  is  due  to  the  fact  that  most  of  the  heat  transferred  through  this  mix 
is  through  the  larger  aggregate  (up  to  38  mm),  and  that  after  the  mineral 
aggregate  voids  are  filled,  the  asphalt  thickness  around  the  larger  particles  is 
not  increased  substantially.  The  continuous  increase  in  k  with  asphalt  content 
may  suggest  a  rearrangement  of  the  mineral  structure  such  that  although  the  mean 
distance  between  particles  has  increased  (as  suggested  by  a  slight  decrease  in 
density),  the  nearest  proximity  points  of  individual  aggregates  may  have  de¬ 
creased. 

The  open  graded  surface  course  mix  with  both  the  limestone  and  the  granite- 
gneiss  aggregate  experienced  first  a  decrease  then  an  increase  in  thermal  con¬ 
ductivity  as  the  asphalt  content  was  increased  from  3.5%  to  6.5%.  This  can  be 
explained  by  recognizing  that  there  are  two  opposing  effects  of  adding  asphalt 
to  a  uniform  aggregate  size.  The  first,  as  mentioned  previously,  is  to  replace 
air  in  the  mineral  aggregate  voids  with  asphalt.  This  would  tend  to  increase 
the  conductivity  of  the  mix--as  seen  in  both  of  the  open  graded  mixes  as  the 
asphalt  content  is  increased  from  5.0%  to  6.5%  (Figure  I - 10 ) .  Secondly,  and 
working  counter  to  the  first  effect,  is  the  influence  of  increasing  the  asphalt 
coating  around  the  particles  as  the  asphalt  content  is  increased.  Since  the 
particles  are  moved  farther  apart  by  the  less  conductive  asphalt,  the  conduc¬ 
tivity  of  the  mix  would  be  reduced--as  seen  in  both  of  the  open  graded  mixes  as 
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the  asphalt  content  is  increased  from  3.5%  to  5.0%  (Figure  I - 10 ) .  The  increase 
in  asphalt  coating  thickness  as  asphalt  content  increases  would  be  greater  in 
the  open  graded  surface  course  mix  than  in  the  other  aggregate  mixes  studied 
because  of  the  smaller  specific  surface  (surface  area  per  unit  volume)  of  the 
open  graded  mix.  The  specific  surface  of  the  open  graded  surface  course  mix  is 
about  40%  of  that  of  the  dense  graded  surface  course  mix. 

The  granite-gneiss  open  graded  mix  had  thermal  conductivity  values  about 
30%  less  than  that  of  the  limestone  open  graded  mix.  This  lower  value  is  the 
result  of  the  lower  conductivity  of  the  granite-gneiss  mineral  (See  Table  1-2). 

The  thermal  conductivity  of  the  dense  graded  surface  course  mix  with  the 
lightweight  aggregate  was  about  40%  less  than  that  of  the  dense  graded  mix  hav¬ 
ing  limestone  aggregate  and  about  11%  less  than  that  of  the  dense  graded  mix 
having  granite-gneiss  aggregate.  This  mix  consisted  (by  weight)  of  25%  expanded 
shale  (from  5  mm  to  13  mm  size),  30%  limestone  aggregate  and  filler,  and  45% 
river  sand  (quartz).  The  expanded  shale  controlled  the  conductivity  of  the  mix 
and  the  asphalt  content  had  little  effect  on  the  conductivity  (Figure  1-10).  It 
is  postulated  that  the  air  trapped  in  the  expanded  shale  is  an  important  limit¬ 
ing  factor  controlling  the  conductivity  of  this  mix. 

Since  the  specific  heat  of  asphalt  is  much  higher  than  that  of  air,  lime¬ 
stone,  granite-gneiss,  and  shale  (Tables  1-1  and  1-2)  it  was  expected  that  as 
the  asphalt  content  of  the  mixes  was  increased  from  3.5%  to  6.5%,  there  would  be 
a  corresponding  increase  in  the  specific  heat  of  the  mixes.  However,  the  spe¬ 
cific  heat  of  the  mixes  is  more  dependent  on  the  aggregate  gradation  and  miner¬ 
alogy  than  the  asphalt  content  (Figure  I - 11 ) .  Based  on  these  data,  it  seems 
that  an  average  value  of  c  for  each  mix,  independent  of  asphalt  content,  could 
be  used. 

The  average  specific  heat  of  the  lightweight  aggregate  mixes  was  much  higher 
than  the  limestone  aggregate  mixes.  This  is  primarily  due  to  the  decrease  in 
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density  of  the  expanded  shale. 

It  is  interesting  to  note  that  the  product  pc  of  the  lightweight  mixes  and 
of  both  dense  and  open  graded  surface  course  limestone  mixes  is  very  nearly  the 
same  (2.5*106  J/m3°C).  This  may  indicate  that  the  specific  heat  is  dependent 
upon  the  density  for  these  mixes.  This  product  for  the  limestone  aggregate  base 
course  mix  is  about  13%  less,  while  the  product  for  both  granite-gneiss  mixes  is 
about  24%  less. 

The  di ffusivity-asphalt  content  relationship  for  each  of  the  mixes  is  shown 
in  Figure  1-12.  For  the  open  graded  surface  course  mixes  (limestone  and  granite- 
gneiss  aggregates),  a  is  essentially  independent  of  asphalt  content;  for  the 
limestone  base  course  mix  and  the  limestone  and  granite-gneiss  dense  graded  sur¬ 
face  course  mixes,  the  diffusivity  varied  about  12%  over  the  range  of  asphalt 
contents  used  in  the  mixes.  The  values  shown  in  Figure  1-12  for  the  limestone 
and  granite-gneiss  mixes  are  well  within  the  range  reported  by  others  (Table 
1-1).  Due  to  the  lower  thermal  conductivity  values,  the  mixes  made  with  the 
granite-gneiss  aggregate  had  diffusivity  values  less  than  that  of  the  comparable 
mixes  made  with  limestone  aggregate.  The  granite-gneiss  dense  graded  surface 
course  mix  had  diffusivity  values  about  9%  less  than  the  limestone  dense  graded 
mix,  while  the  diffusivity  of  the  granite-gneiss  open  graded  mix  was  about  25% 
less  than  that  of  the  limestone  open  graded  mix.  The  diffusivity  of  the  dense 
graded  surface  course  using  lightweight  aggregate  was  much  less  than  that  of  the 
dense  graded  limestone  mix  and  was  found  to  be  independent  of  asphalt  content. 
Since  the  product  pc  for  the  lightweight  mixes  is  similar  to  that  of  the  lime¬ 
stone  mixes,  the  reduction  in  a  is  due  to  the  smaller  thermal  conductivity  of 
the  lightweight  mixes  (Figure  I - 10 ) . 

CHAPTER  I-F  RESULTS  OF  PROPERTY  ANALYSIS  OF  FIELD  SPECIMENS 

This  chapter  presents  the  physical  and  thermal  analyses  of  the  asphaltic 
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concrete  pavement  cores  taken  from  Myrtle  Beach  and  Loring  Air  Force  Bases. 
These  bases  were  selected  due  to  their  geographical  extremes  (South  Carolina  and 
northern  Maine)  which  would  insure  different  asphaltic  concrete  pavement  mix 
designs. 

Seven  10.2  cm  diameter  cores  were  drilled  at  different  locations  on  both 
the  Myrtle  Beach  and  Loring  Bases.  The  top  six  centimeters  of  each  core  was 
sliced  with  a  rock  saw  to  make  three  disks  about  two  centimeters  thick  each. 
The  density  of  each  disk  was  determined  according  to  ASTM  D1188,  and  thermo¬ 
couple  recessing  grooves  were  sawn  into  one  surface  of  each  disk.  The  thermal 
properties  of  each  disk  were  then  obtained  by  the  procedure  described  previously 
(Chapter  I-C). 

One  of  the  three  disks  from  each  core  was  used  in  an  asphalt  extraction 
test,  ASTM  D2172-81,  to  determine  the  amount  of  asphalt  cement  present  in  the 
different  pavements.  Also  a  grain  size  analysis  was  conducted  on  the  aggregate 
remaining  from  each  of  the  extraction  tests. 

A  summary  of  the  result  of  the  thermal  and  physical  analyses  of  the  field 
pavement  samples  is  shown  in  Figures  1-13  and  1-14  and  Tables  1-5,  1-6,  and  1-7. 

In  the  analysis  of  the  data,  the  mean  and  standard  deviation  of  the  density  (p), 
percent  asphalt  by  weight,  diffusivity  (a),  conductivity  (k),  and  the  specific 
heat  (c)  are  calculated  to  indicate  the  average  values  and  the  variations,  respec¬ 
tively  (Table  1-5).  Tables  1-6  and  1-7  show  the  thermal  and  physical  properties 
of  the  three  groups  of  similar  pavements  from  Myrtle  Beach  and  Loring  Air  Force 
Bases.  A  group  is  defined  as  samples  from  locations  on  the  base  with  similar 
histories,  features,  and  pavement  thicknesses.  The  features  and  history  of 
Loring  Group  A  are  very  well  defined  since  all  the  pavements  in  this  group  were 
resurfaced  in  the  summer  of  1982.  The  Myrtle  Beach  Group  and  Loring  Group  B 
each  have  similar  histories  but  the  features  and  pavement  thicknesses  are  not  as 
closely  matched  as  Loring  Group  A. 
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rigure  1-13.  Range  of  Grain  S i /e  Distributions  for  the  Luring  Samples 


Table  1-5.  Summary  of  Physical  and  Thermal  Properties  of  Myrtle  Beach 
(M)  and  Loring  (L)  Air  Force  Base  Field  Samples. 


ID# 

(Mg/m3) 

%  Asphalt 

PROPERTIES 

(average  of  3  samples) 

k 

(w/m°C) 

c 

( J/kg°C) 

a 

(m2/s*10-7 ) 

Ml 

2.26 

4.8 

0.80 

425 

8.4 

M2 

2.25 

4.6 

0.71 

412 

7.6 

M3 

2.22 

8.0 

0.93 

500 

6.3 

M4 

2.28 

8.0 

0.82 

561 

6.4 

M5 

2.27 

6.  7 

0.93 

572 

7.2 

M6 

2.22 

6.7 

0.63 

445 

6.5 

M7 

2.37 

6.7 

0.67 

423 

6.9 

Mean 

2.27 

6.5 

0.78 

477 

7.0 

S.D. 

0.05 

1.4 

0. 12 

68 

0.8 

LI 

2.41 

8.9 

0.68 

501 

5.5 

L2 

2.38 

8.9 

0.80 

636 

5.3 

L3 

2.32 

8.9 

0.83 

651 

5.6 

L6 

2.38 

8.5 

0.81 

629 

5.4 

L7 

2.49 

5.7 

0.90 

763 

5.0 

L8 

2.45 

7.0 

0.80 

556 

5.9 

L9 

2.42 

7.0 

0.82 

595 

5.7 

Mean 

2.41 

7.8 

0.81 

618 

5.5 

S.D. 

0.05 

1.3 

0.07 

82 

0.3 
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Table  1-6.  Summary  of  Physical  and  Thermal  Properties  of 
Myrtle  Beach  Group. 


PROPERTIES 

(average  of  3  samples) 


ID#  p  %  Asphalt  k  c  a 

(Mg/m3)  (w/m°C)  (J/kg°C)  (m2/s*10-7) 


M3 

2.22 

8.0 

0.93 

500 

6.3 

M4 

2.28 

8.0 

0.82 

561 

6.4 

M7 

2  37 

6.7 

0.67 

423 

6.9 

Mean 

2.29 

7.6 

0.81 

495 

6.5 

S.D. 

0.08 

0.7 

0.13 

69 

0.3 

Table  1-7.  Summary  of  Physical  and  Thermal  Properties  of 
Loring  Group  A  and  Group  B. 


ID# 

(Mg/m3) 

%  Asphalt 

PROPERTIES 

(average  of  3  samples) 

k 

(w/m°C) 

c 

(J/kg°C) 

a 

(m2/s*10  7) 

Group  A 

LI 

2.41 

8.9 

0.68 

501 

5.5 

L2 

2.38 

8.9 

0.80 

636 

5.3 

L3 

2.32 

8.9 

0.83 

651 

5.6 

L4 

2.38 

8.5 

0.81 

629 

5.4 

Mean 

2.37 

8.8 

0.78 

604 

5.5 

S.D. 

0.04 

0.2 

0.07 

69 

0.1 

Group  B 

L7 

2.49 

5.7 

0.90 

763 

5.0 

L8 

2.45 

7.0 

0.80 

556 

5.9 

L9 

2.42 

7.0 

0.82 

595 

5.7 

Mean 

2.45 

6.6 

0.84 

638 

5.5 

S.D. 

0.04 

0.7 

0.05 

110 

0.5 
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The  extraction  test  results  indicate  a  range  of  asphalt  contents  from  4.6% 
to  8.0%  (Table  1-5).  The  Myrtle  Beach  samples  have  an  average  asphalt  content 
of  about  6.5%,  while  the  Loring  samples  have  a  higher  average  of  7.8%.  The 
course  aggregate  in  the  Myrtle  Beach  samples  is  granite,  and  the  fine  aggregate 
is  quartz  sand.  The  Loring  pavement  aggregate  is  composed  of  a  mixture  of  slate 
and  basalt  (10)  as  both  the  coarse  and  fine  aggregate. 

Grain  size  distributions  are  shown  in  Figures  1-13  and  1-14.  The  dashed 
line  on  each  of  the  figures  is  the  grain  size  distribution  curve  for  the  dense 
graded  surface  course  mix  which  was  presented  in  Chapter  I-D.  The  range  of  grain 
sizes  from  the  Loring  samples  very  closely  resembles  the  laboratory  dense  graded 
surface  course  curve  (Figure  1-13).  The  Myrtle  Beach  samples  contain  a  higher 
percentage  of  fine  aggregate  (sand  size)  than  the  laboratory  dense  graded  sur¬ 
face  course  mix  (Figure  1-14). 

The  density  of  each  mix  is  shown  in  Table  1-5.  The  samples  from  Myrtle 
Beach  have  densities  averaging  about  6%  less  than  the  samples  from  Loring.  The 
mean  density  of  each  of  the  three  groups  (Tables  1-6  and  1-7)  is  within  2%  of 
the  mean  of  their  parent  data  sets  (Table  1-5). 

The  thermal  conductivity  of  each  mix  is  shown  in  Table  1-5.  Each  value 
represents  the  average  of  three  samples.  The  mean  value  of  each  data  set  is 
very  nearly  the  same,  with  the  Loring  set  having  a  slightly  higher  average  of 
thermal  conductivity.  This  is  probably  due  to  the  lower  conductivity  of  the 
quartz  sand  fine  aggregate  present  in  the  Myrtle  Beach  pavements  (see  Table 
1-2).  There  is  less  variation  in  the  Loring  samples  as  shown  by  the  smaller 
standard  deviation  value.  The  Myrtle  Beach  Group  (Table  1-6)  has  a  mean  thermal 
conductivity  value  about  4%  higher  that  the  mean  conductivity  of  all  the  Myrtle 
Beach  data  (Table  1-5).  The  Loring  Groups  have  mean  conductivity  values  of  4% 
less  (Group  A)  and  4%  greater  (Group  B)  than  the  mean  conductivity  of  all  the 
Loring  data  (see  Tables  1-5  and  1-7). 


The  average  specific  heat  values  (shown  in  Table  1-5)  for  the  Myrtle  Beach 
samples  are  much  less  than  the  average  values  of  the  Loring  samples  (about  22% 
less).  This  is  due  to  the  higher  specific  heat  of  the  Loring  slate-basalt 
aggregate  than  that  of  the  granite  and  quartz  aggregate  in  the  Myrtle  Beach 
pavements  (See  Table  1-2).  The  mean  specific  heat  values  of  each  of  the  three 
groups  (Tables  1-6  and  1-7)  are  within  4%  of  the  mean  specific  heat  value  of 
their  parent  data  sets  (Table  1-5). 

The  average  thermal  diffusivity  (Table  1-5)  of  the  Loring  samples  is  much 
(21%)  less  than  that  of  the  Myrtle  Beach  samples  since  both  the  density  and  the 
specific  heat  of  the  Loring  samples  have  higher  values.  The  mean  diffusivity  of 
both  of  the  Loring  Groups  (Table  1-7)  is  the  same  as  that  of  the  entire  Loring 
data  set,  while  the  mean  diffusivity  of  the  Myrtle  Beach  Group  (Table  1-6)  is  7% 
less  than  that  of  all  the  Myrtle  Beach  data. 

Since  the  uncertainty  of  each  of  the  measured  thermal  properties  is  about 
5.5%  (see  Appendix  I-B),  it  can  be  said  that  each  group's  mean  thermal  property 
values  (Tables  1-6  and  1-7)  is  essentially  the  same  as  its  parent  data  set 
(Table  1-5),  i.e.  the  uncertainty  of  the  measured  values  is  greater  than  the 
variation  of  the  mean  values.  Also,  the  mean  conductivity  of  each  of  the  Loring 
and  Myrtle  Beach  data  sets  (Table  1-5)  is  essentially  the  same,  since  the  per¬ 
cent  difference  between  the  two  mean  values  is  less  than  the  percent  uncertainty 
(Appendix  I-B).  However,  the  mean  specific  heat  and  diffusivity  of  the  Loring 
and  Myrtle  Beach  data  sets  are  quite  different  due  to  the  different  densities 
and  aggregate  types. 
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PART  II 

A  NUMERICAL  MODEL  OF  THE 
HEATING  PROCESS  FOR 
ASPHALTIC  CONCRETE 


PAVEMENT 


CHAPTER  1 1 -A  INTRODUCTION 


The  process  of  heating  an  asphalt  pavement  prior  to  its  removal  has  several 
inherent  problems  that  have  not  yet  been  adequately  resolved.  Foremost  among 
these  problems  is  the  inadvertent  depletion  of  plasticizing  agents  in  the  as¬ 
phalt.  These  volatile  components  are  driven  off  when  the  temperature  of  the 
asphalt  exceeds  a  threshold  value  of  approximately  505  (°K)  during  the  heating 
process.  This  loss  of  plasticizing  agents  is  harmful  since  it  is  known  to  pro¬ 
mote  early  aging  of  asphaltic  pavements.  Since  the  heating  is  usually  accom¬ 
plished  with  fossil  fuel-fired  heaters,  the  chemical  reaction  of  certain  com¬ 
bustion  products  with  the  asphalt  also  presents  a  problem.  These  and  other 
problems  associated  with  the  heating  of  asphalt  pavements  have  motivated  the 
construction  of  the  analytical  model  of  the  heating  process  presented  in  this 
report. 

CHAPTER  II-B.  A  TYPICAL  STRATEGY  FOR  HEATING  ASPHALTIC  PAVEMENT  PRIOR  TO 
REMOVAL  OF  A  SURFACE  LAYER 

As  stated  above,  the  pavement  must  be  heated  to  its  removal  depth  in  order 
to  reduce  its  shearing  strength.  Specifically,  the  object  of  this  heating  pro¬ 
cess  is  to  raise  the  temperature  of  the  entire  layer  of  pavement  to  be  removed 
to  above  a  minimum  desirable  temperature  for  removal  =  363°K)  without 
exceeding  the  maximum  allowable  temperature  of  505°K  anywhere  in  the  pavement. 

For  economic  reasons,  the  heating  process  must  be  completed  as  rapidly  as 
possible.  Due  to  the  extremely  low  thermal  diffusivity  of  asphaltic  pavement, 
however,  a  typical  pavement  heater  is  capable  of  raising  the  surface  temperature 
of  the  pavement  to  above  the  maximum  allowable  value  (505°K)  long  before  the 
temperature  at  the  removal  depth  has  reached  the  minimum  desirable  temperature 
for  removal  (363°K).  Thus,  some  paving  contractors  have  adopted  a  simple  strat- 


egy  for  heating  asphaltic  pavement  that  may  be  summarized  in  the  following  steps: 

1)  The  pavement  is  heated  until  the  maximum  allowable  temperature  is  reached 
at  the  upper  surface. 

2)  Heating  is  then  discontinued  and  the  surface  is  insulated  for  an  arbitrar¬ 
ily  specified  period  of  time  during  which  the  surface  temperature  decreases 
due  to  the  diffusion  of  thermal  energy  into  the  interior  of  the  pavement  by 
conduction. 

3)  The  insulation  is  removed  from  the  surface  and  the  heating  is  begun  again. 
This  heating-insulation  cycle  is  repeated  until  the  temperature  reaches  the 
desired  value  (363°K)  at  the  depth  which  is  to  be  removed. 

The  analytical  model  of  the  pavement  heating  process  presented  in  this 
report  is  based  on  the  heating  strategy  described  above. 

CHAPTER  II-C.  A  DETAILED  DESCRIPTION  OF  ASPHALTIC  PAVEMENT 
An  asphaltic  pavement  is  composed  of  solid,  semisolid,  or  gaseous  phases. 
The  solid,  or  aggregate,  phase  typically  consists  of  sand,  gravel,  crushed  stone, 
slag  and  mineral  filler.  The  semi-solid  phase  is  a  viscoelastic  asphaltic  mate¬ 
rial  which  is  derived  from  petroleum  and  which  may  have  widely  varying  proper¬ 
ties.  Approximately  5%  of  the  mixture  by  weight  is  composed  of  the  asphaltic 
material  which  serves  to  bind  together  the  aggregate.  The  gaseous  phase  is  com¬ 
posed  of  the  air  which  fills  the  voids.  Water  in  either  vapor  or  liquid  form 
may  also  occupy  these  voids. 

A  typical  asphaltic  pavement  is  constructed  of  three  layers,  or  "courses", 
as  shown  in  Figure  II-l.  The  uppermost,  or  "wearing",  course  is  a  well-compact¬ 
ed,  hot-rolled  asphaltic  mix  which  has  a  very  low  permeability  to  water,  while 
the  lowest,  or  "base",  course  may  have  a  high  percentage  of  voids  that  are  eas¬ 
ily  penetrated  by  water.  The  middle,  or  "binder,"  course  is  usually  of  an  inter¬ 
mediate  composition.  Measurements  of  the  thermophysical  properties  of  asphaltic 
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concrete  have  been  reported  by  Kavianipour  and  Beck  (19)  and  in  Part  I  of  this 
report. 

CHAPTER  II-D.  ANALYSIS  OF  THE  HEATING  PROCESS 

Problem  Statement 

It  is  desired  to  develop  an  analytical  model  of  the  transient  heat  transfer 
that  occurs  during  the  heating  of  an  asphaltic  pavement  by  the  combined  action 
of  convection  and  radiation  at  the  upper  surface  of  the  wearing  course. 

Modeling  the  Composite  Pavement  as  a  Homogeneous  Body 

The  heat  transfer  within  the  pavement  is  entirely  by  conduction.  Since  the 
pavement  is  composed  of  three  courses,  each  with  different  thermophysical  pro¬ 
perties,  it  would  seem  that  the  pavement  should  be  modeled  as  a  three- layer  com¬ 
posite  body.  The  variations  between  the  properties  of  the  three  courses  are 
relatively  small,  however,  which  suggests  that  the  pavement  might  be  adequately 
described  by  modeling  it  as  a  homogeneous  body.  Thus,  an  analytical  study  was 
performed  to  justify  the  use  of  this  approximation.  This  study,  which  is  des¬ 
cribed  in  detail  in  Appendix  I I- A ,  shows  that  the  errors  incurred  by  modeling  a 
three-layer  asphaltic  pavement  as  a  homogeneous  (one-layer)  body  are  negligibly 
small.  Therefore,  the  analytical  model  of  the  asphaltic  pavement  heating  pro¬ 
cess  presented  in  this  report  is  based  on  the  assumption  that  an  actual  three- 
layer  pavement  may  be  accurately  modeled  as  a  homogeneous  body. 

The  Analytical  Model 

A  schematic  of  the  idealized  heater-pavement  system  being  analyzed  is  shown 
in  Figure  1 1-2.  The  analytical  model  of  this  system  is  based  on  the  following 
assumptions: 

(1)  The  three-course  composite  pavement  may  be  accurately  modeled  as  a  homo¬ 
geneous  (one- layer)  body  with  constant  thermal  properties. 
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(2)  The  heat  transfer  within  the  pavement  takes  place  entirely  by  transient, 
one-dimensional  conduction. 

(3)  The  heat  transfers  along  the  upper  surface  of  the  wearing  course,  the  com- 
busion  gases,  and  the  bottom  surface  of  the  heater  occur  by  combined  radia¬ 
tion  and  convection. 

Note:  The  role  of  convection  in  the  overall  heat  transfer  process  at  the 
upper  surface  of  the  pavement  is  currently  not  well-understood.  Several 

models  depicting  the  possible  contributions  of  convection  to  this  heat 

transfer  process  are  briefly  discussed  in  Appendix  II-B. 

(4)  The  radiative  exchange  between  the  upper  surface  of  the  pavement  and  the 

lower  surface  of  the  heater  is  the  same  as  the  radiative  exchange  between 

two  infinite,  parallel  gray  planes.  (The  gap  between  these  two  planes  is 
assumed  to  be  filled  with  a  nonparticipating  medium  since  the  degree  of 
participation  of  the  hot  combustion  gases  in  the  radiative  exchange  is  also 
presently  not  well  understood.) 

(5)  The  temperature  at  the  bottom  surface  of  the  pavement  remains  constant  and 
is  equal  to  the  soil  temperature. 

(6)  The  initial  temperature  distribution  in  the  pavement  is  known  and  the  heat¬ 
ing  process  is  initiated  by  exposing  the  upper  surface  of  the  pavement  to 
the  lower  surface  of  the  heater  (as  shown  in  Figure  1 1 -2 ) . 

The  above  assumptions  yield  the  following  governing  equation  for  the  tran¬ 
sient  conduction  in  the  pavement  along  with  its  initial  and  boundary  condi¬ 
tions: 

3T*  =  aV  (0  <  x*  <  L)  (la) 

at*  ax*2 

T*(x*,0)  =  T*(x*)  (0  <  x*  <  L)  (lb) 
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T*(L,t*)  =  T*oi-j  (=  constant) 


[T*4-  T*4(0,t*)]  +  hc[T*  -  T*(o,t*)] 


(lc) 


=  -k  OT* 
9x* 


(Id) 


where: 

T*  s  pavement  temperature 

x*  =  spatial  coordinate  in  the  pavement 

t*  =  time 

L  =  overall  thickness  of  the  pavement 

T*oii=  temperature  at  the  bottom  surface  of  the  pavement 


T* 

H 


heater  surface  temperature 


T* 

R 


a 

k 

e_ 


an  appropriate  reference  temperature  for  the  convective  heat 
transfer 

convective  heat  transfer  coefficient 
thermal  diffusivity  of  the  pavement 
thermal  conductivity  of  the  pavement 
emissivity  of  the  upper  surface  of  the  pavement 
emissivity  of  the  heater  surface 

_  O 

Stefan  Boltzman  constant  =  (5.6697  *  10  )  /  w  \ 


The  transient  conduction  problem  defined  by  equations  (la)-(ld)  is  non 
linear  because  of  the  radiation  boundary  condition  at  the  upper  surface.  The 


use  of  a  so-called  "radiation  heat  transfer  coefficient"  to  linearize  this  bound¬ 
ary  condition  is  inappropriate  in  this  problem  because  of  the  large  temperature 
differences  between  the  pavement  and  heater  surfaces.  Thus,  the  problem  will  be 
solved  numerically  using  finite  difference  techniques. 

Formulation  of  the  Numerical  Model 

For  the  purposes  of  the  numerical  model,  the  pavement  is  represented  by  (M 
+  1)  nodes  as  shown  in  Figure  1 1- 3 .  The  two  surface  nodes  (nodes  1  and  M  +  1) 
each  represent  a  layer  of  pavement  with  a  thickness  of  Ax*/2  (where  Ax*  =  L/M), 
while  the  (M  -  1)  interior  nodes  each  represent  a  layer  of  pavement  of  thickness 
Ax*.  Since  the  temperature  at  the  bottom  of  the  slab  is  assumed  to  be  a  given 
constant  (T*  . the  temperature  of  the  (M  +  l)-th  node  is  known.  Thus,  only 
M  equations  must  be  written  to  determine  the  temperatures  of  the  remaining  nodes 
(nodes  1  through  M). 

An  energy  balance  on  the  i-th  interior  node  gives 


''Rate  of  conduction^ 
heat  transfer  to 
vnode  i  / 


;Rate  of  conduction' 
heat  transfer 
^from  node  i  > 


^Rate  of  change  of  v 
thermal  energy  stored] 
<in  node  i  / 


r T?-x  -  T?  i 

rT*  - 

-  -  IrA  1 

T*  *  i] 

Lx*  ->•*  J 

—  KM  1 

Lx*  - 

X  *  +  1  J 

which  reduces  to 


dT* 

—L~  ~  - -  [T*  -  2T?  +  T*  .]  (2  <  i  <  M)  (la) 

dt*  (Ax*;  11  1  11 

where: 


a  =  k_  =  thermal  diffusivity  of  the  pavement 
pc 
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heater  surface 


node  1 


T  *(  0  ,  t  J 


\  ,r  „ 


X*  -o 


•  i -  1  th  node 


• l  th  node 


■i+i  th  node 


genera ] 
interior 


/  7  / 

^  ^  M+l  th  node  /  /  ^ 

Figure  II  -3. 

Nodal  Representation  Of  The  Pavement  For  The 
Finite-Difference  Model 

Similarly,  an  energy  balance  for  node  1  (the  node  at  the  upper  surface  of 
the  pavement)  while  the  pavement  is  being  heated  gives: 


(Rate  of  radiation  heat 
transfer  to  node  1 


/Rate  of  convection  heat 
(transfer  to  node  1  t 


(Rate  of  radiation\ 
heat  transfer  I  + 
from  node  1  / 


(Rate  of  conduction 
heat  transfer 
from  node  1 


)/ Rate  of  change  of\ 
+  (  thermal  energy  ) 

\stored  in  node  1  / 


or, 


°  »h-i  (th4  -  t*4)  *  hc  (tr  -  t;> 


+  pcAAx* 
2 


dT^ 

dF" 


which  reduces  to 


[a  yH.1  (T*4  -  T*4)  *  hc  (T«  -  T*) 


LAJ 


Ax’ 


(T*  -  T*)] 


(lb) 


where: 

Yh-1  2  radiation  heat  transfer  modulus  that  accounts  for  the  emissivities 
and  relative  geometries  of  the  pavement  and  heater  surfaces 

(=  1/e  ♦  l/£u  -  1  } 

s  n 

(x*  -  xp  =  -  (Ax* ) 

Note:  The  numerical  model  must  be  capable  of  computing  the  transient 
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temperature  distributions  in  the  pavement  during  both  the  heating  and  in¬ 
sulation  periods  of  the  cyclic  process  described  earlier  in  this  report. 


Thus,  an  energy  balance  on  node  1  while  the  surface  is  insulated  gives 


"Rate 

of  conduction 

Rate  of  change  of 

heat 

transfer  from 

+ 

thermal  energy 

=  0 

.node 

1 

stored  in  node  1  . 

=  0 


which  reduces  to: 


dT* 

dt* 


2a 


(A  x*Y 


[T*  -  T*] 


(lc) 


The  initial  condition  is  given  by  specifying  the  initial  temperature 


for  [1  <  i  <  (M  +  1)]  at  t*  =  0 


(Id) 


at  each  nodal  point  (where  C*  is  an  appropriate  set  of  constants).  Finally,  the 
boundary  condition  at  the  bottom  surface  of  the  pavement  is  given  by 

TS.l  *  T;„f,  (for  t*  >  »)•  <le> 


In  summary,  the  numerical  model  (written  in  terms  of  physical  variables) 
consists  of  equations  (la)  through  (le).  These  equations  contain  a  number  of 

independent  parameters  (i.e.,  a,  h  .  T£,  e  ,  eu,  . ).  As  in  all  heat  trans- 

fer  problems,  it  is  convenient  to  reduce  the  number  of  independent  parameters 
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enhances  our  physical  understanding  of  the  solution,  but  also  greatly  reduces 
the  number  of  computations  required  to  complete  a  parametric  study  of  the  pro¬ 
blem  by  presenting  the  solutions  in  their  most  general  forms.  Thus,  the  follow¬ 
ing  set  of  nondimensional  variables  are  introduced 


T  s  T* 

T*~ 

1  nr 


(L2/a) 


where: 


the  minimum  temperature  to  which  the  pavement  must  be  heated  at 
the  depth  to  which  it  is  to  be  removed  before  the  removal 
can  take  place  (=  363°K) 


Note: 


(1)  The  temperature  must  be  nondimensional ized  as  shown  above  instead  of  in 
the  usual  manner;  that  is,  by  letting 

T  -  T„rr 


because  of  the  nonlinear  radiation  boundary  condition. 


(2)  Tjij^  was  chosen  as  the  temperature  scale  because  it  is  the  only  tem¬ 
perature  in  the  entire  problem  with  a  constant  value. 

Substituting  the  nondimensional  variables  in  equations  (la)  through  (le) 
yields  the  following  set  of  equations  which  comprise  the  nondimensional  form  of 
the  numerical  model. 


CTi-i  -  2  Ti  +  Tm]  (2  <  1  <  M> 
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—  =  — [(A)  (Tj  -  T*)  +  (Bi  )(T  -  T  )  -  (  -^  )  (T  -T  )]  (2b) 

dt  Ax  H  1  R  1  Ax  1  2 


—  =  -  [  — -  ]  CT-,  -  T2J 

dl  (AX)2 


T.  =  C.  for  [1  <  i  <  (M  +  1)]  at  t  =  0 


Tm+1  =  Tsoil  (for  1  >  0) 


where: 


oL*h-i  Tem 


aL  J* 


eH  £s 


=  radiation-conduction  parameter 


hcL 

Bi  =  -  =  Biot  number 

k 


=  dimensionless  heater  temperature 


=  dimensionless  reference  temperature  for  convection 


Implementation  of  the  Numerical  Model 

The  nondimensional  form  of  the  numerical  model  given  by  equations  (2a)-(2e) 
was  programmed  for  execution  on  a  digital  computer.  This  program  incorporates 
the  heating  technique  described  earlier  in  this  report;  that  is,  the  heating  is 
begun  and  continues  until  either  the  surface  temperature  exceeds  a  maximum  allow¬ 
able  value 
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or  the  temperature  at 


X* 

A  P£|V| 

XREM  =  -  reaches  the  value  required  for 

removal  [TREM  =  1],  at  which  time  the  process  is  terminated.  If  the  temperature 
at  xREM  has  not  reached  the  desired  value  of  TREM  =  1  when  reaches  (T1)MAX, 
the  heating  is  terminated  and  the  upper  surface  is  insulated  for  an  arbitrarily 
specified  period  of  time, 

t* 

*  _  lINSUL 

^INSUL  = 

(L  /a) 

This  cycle  is  repeated  unti  1  T  =  trem  =  1  at  xREM-  Thus,  specifying  the  values 
of  the  eight  dimensionless  parameters  A,  Bi ,  TH ,  TR1  TREM,  xR£M,  (Tj)^  and 

^INSUL’  a^ong  w^tfl  initial  temperature  distribution,  completely  specifies 
the  problem. 

The  solution  of  the  numerical  model  employs  a  simple,  explict  finite  dif¬ 
ference  technique.  A  brief  description  of  this  technique  and  the  reasons  for 
its  selection  are  given  in  Appendix  II-C.  A  listing  of  the  resulting  FORTRAN 
code,  HEAT1,  is  given  in  Appendix  II-D. 

CHAPTER  II-E.  PRESENTATION  OF  RESULTS  FOR  THE  "NOMINAL"  HEATING  PROCESS  IN  THE 
PRESENT  INVESTIGATION 

To  illustrate  some  typical  output  from  the  HEAT1  code,  representative 

temperature  histories  and  distributions  for  the  "nominal"  heating  case  in  this 

investigation  are  presented  here  in  graphical  form.  The  nominal  heating  case  is 

defined  by  the  values  of  the  eight  dimensionless  parameters  shown  in  Table 
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II-l.  Note  that  this  case  assumes  the  surface  of  the  pavement  is  heated  by 
radiation  alone  and  that  the  initial  temperature  distribu¬ 
tion  in  the  pavement  is  uniform. 

Figure  1 1-4  shows  the  surface  and  removal  depth  temperature  histories  for 
the  nominal  heating  case.  The  monontonical ly  rising  surface  temperature  clearly 
shows  that  no  insulation  periods  are  required  for  the  nominal  case.  The  heater 
temperature  is  low  enough  that  thermal  energy  can  diffuse  into  the  pavement  by 
conduction  fast  enough  to  raise  the  temperature  at  the  removal  depth  to  its  re¬ 
quired  value  (=  1.0)  before  the  surface  temperature  exceeds  the  maximum  allow¬ 
able  value  (=  1.394)  in  the  pavement.  The  Fourier  number  at  the  end  of  the  heat- 

-3 

ing  period  is  (13.4  *  10  ),  which  corresponds  to  a  heating  time  of  31.5  min. 

Figure  II-5  shows  the  dimensionless  temperature  profiles  in  the  pavement  at 
different  dimensionless  times  during  the  heating  period.  These  profiles  dra¬ 
matically  illustrate  the  exceptionally  low  value  of  the  thermal  diffusivity  of 
asphalt.  At  the  end  of  the  heating  period  (31.5  min),  thermal  energy  from  the 
surface  heating  has  diffused  through  only  about  40%  of  the  total  thickness  of 
the  pavement. 
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T.  =  C.  =  T  ..  =  SQl1  =  0.807  For  [1  <  i  <  (M  +  1)]. 

T  1  SOI  I  T* 

1  REM 

Note:  The  nominal  heating  case  assumes  a  uniform  initial  temperature  distribu¬ 
tion. 

Table  II-l.  Values  of  the  Nondimensional  Parameters  in  the  Heatl  Code  for  the 
Nominal  Heating  Case 
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PART  III 

VALIDATION  OF  THE 
THERMAL  MODEL 


CHAPTER  1 1 1 -A  INTRODUCTION 


To  acquire  the  data  needed  validate  the  thermal  model  of  asphaltic  pavement, 
field  tests  were  conducted  at  two  U.  S.  Air  Force  Bases  --  Myrtle  Beach  Air 
Force  Base  (AFB)  in  South  Carolina  and  Loring  AFB  in  northern  Maine.  This  part 
of  the  research  effort  provides  a  description  of  the  heater  tests  and  comments 
on  the  validity  of  the  thermal  model.  The  field  tests  involved  heating  the 
asphaltic  concrete  pavement  using  a  propane  fired  field  heater  and  measuring 
temperatures  at  various  depths  within  the  pavement  as  a  function  of  time.  To 
determine  the  validity  of  the  model  presented  in  Part  II,  the  thermal  proper¬ 
ties  of  the  appropriate  field  sample  (See  Chapter  I  -  F ;  emissivity  values  for 
steel  and  asphalt  were  assumed)  were  used  in  the  numerical  model  of  the  heating 
process  along  with  pertinent  field  heater  characteristics  to  predict  the  tem¬ 
poral-spatial  temperature  history  of  an  in-service  asphaltic  concrete  pavement 
at  the  site  the  sample  was  obtained.  The  predicted  temperature-time  history  at 
various  depths  in  the  pavement  was  then  compared  to  the  temperatures  measured  at 
a  particular  site  during  a  field  test.  Good  agreement  between  observed  (mea¬ 
sured)  and  predicted  temperature  data  would  establish  the  validity  of  the  model. 
This  part  of  the  research  effort  provides  a  description  of  the  heater  tests  and 
comments  on  the  validity  of  the  thermal  model. 

CHAPTER  I I 1-8  DESCRIPTION  OF  FIELD  HEATER  TESTS  AND  DISCUSSION  OF  RESULTS 

Seven  sites  were  selected  at  each  Base,  and  15.2  cm  cores  were  drilled 
vertically  through  the  pavement  at  each  site.  From  inside  the  core  hole,  four 
lateral  holes  (0.32  cm  in  diameter)  and  about  5  cm  long  were  drilled  horizon¬ 
tally  into  the  pavement  which  provided  thermocouple  placement  at  depths  within 
the  pavement  varying  from  about  0.5  cm  to  about  5.6  cm  (See  Figure  1 1 1  - 1 ) .  The 
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Figure  III-l.  Plan  and  Section  View  of  a  Typical  Test  Hole 


lateral  holes  were  filled  with  high  conductivity  silicon  grease  (Omegatherm  201 
from  Omega  Engineering,  Inc.),  and  a  24  AWG  cooper-constantan  thermocouple  was 
inserted  into  each  lateral  hole.  At  the  pavement  surface,  a  24  AWG  copper-con- 
stantan  thermocouple  was  embedded  in  a  globule  of  high  conductivity  silicon 
grease.  These  five  thermocouples  (and  a  sixth  which  was  attached  to  the  heater) 
led  to  the  same  temperature  measuring  system  (Apple  11+  computer  and  TEMPSENSE 
software)  used  in  the  thermal  property  analysis  described  in  Chapter  I-C. 

A  schematic  of  the  field  heater  is  shown  in  Figure  I I 1-2 .  The  heater  was 
manufactured  by  Aeroil  Products  Co.  ,  Inc.  of  New  Jersey  and  is  designed  to  be 
used  in  pavement  pothole  patching.  Although  it  was  understood  that  with  this 
heater  the  surface  heating  of  asphaltic  pavement  could  be  modeled  as  a  purely 
radiation  process  (26),  test  data  indicated  that  convective  cooling  around  the 
edges  of  the  heater  was  an  important  characteristic  and  therefore  was  incorpo¬ 
rated  in  the  model.  A  20  AWG  chromel-alumel  heater  thermocouple  was  mounted 
between  the  two  center  heater  tubes  at  the  same  height  as  the  bottom  of  the 
heater  tubes  (See  Figure  I I I -1) .  It  was  found  later  that  this  location  gave 
temperature  readings  that  were  not  representative  of  the  temperature  of  the 
burner  tubes.  Heater  temperatures  were  corrected  as  described  below. 

Heater  tests  were  carried  out  with  thermocouples  attached  to  different 
points  on  the  heater.  Since  the  temperature  of  the  burner  tubes  was  important 
in  the  analysis,  a  correlation  was  established  between  the  burner  tube  tempera¬ 
ture  and  the  temperature  at  a  point  between  the  center  heater  tubes.  The  field 
neater  temperature  data  were  then  corrected  using  the  adjusted  temperature. 
This  adjustment  was  critical  in  the  analysis  because  the  radiant  heat  delivered 
to  the  pavement  surface  is  a  function  of  the  fourth  power  of  the  heater  tempera- 
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ture  (Tj_j  )  and  small  errors  in  T^  would  lead  to  very  large  errors  in  the  ana¬ 
lysis. 

The  Myrtle  Beach  field  tests  were  run  in  March,  1983.  The  coastal  winds 
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Figure  1 1 1  -  2 .  Plan,  Profile,  and  Section  View  of  Field  Heater. 


prevailing  at  the  time  disturbed  the  temperature  patterns  of  the  air  between 
the  heater  and  the  pavement  surface  and  all  the  heat  generated  by  the  heater  was 
not  delivered  to  the  pavement  surface  beneath  the  heater.  Thus,  for  a  given 
heater  temperature,  the  temperature  of  the  pavement  surface  fluctuated  widely 
with  the  wind  velocity  (both  speed  and  direction).  Hence,  the  temperature  data 
were  considered  to  be  unreliable  and.no  meaningful  analysis  could  be  carried  out 
on  the  Myrtle  Beach  data. 

The  problem  was  solved  at  Loring  Air  Force  Base  when  these  tests  were  run 
in  June,  1983  by  erecting  a  corrugated  steel  fence  (wind  fence)  about  3  m  square 
and  about  0.8  m  high  around  the  heater  and  the  test  hole.  That  the  wind  fence 
was  very  effective  in  reducing  the  wind  velocity  between  the  heater  and  the  pave¬ 
ment  surface  was  shown  by  the  less  variable  pavement  surface  and  heater  tempera¬ 
ture  data. 

A  typical  test  proceeded  as  follows: 

o  A  15.2  cm  diameter  core  and  lateral  thermocouple  holes  were  drilled, 
o  The  thermocouple  hole  depths  were  measured  and  recorded, 
o  The  wind  fence  was  set  up  to  enclose  the  test  hole  and  the  heater, 
o  An  asbestos  board  was  placed  underneath  the  heater  to  prevent  scor¬ 
ing  the  pavement  while  the  heater  was  reaching  steady  state, 
o  The  thermocouple  holes  were  filled  with  high  conductivity  silicon 
grease  and  the  four  subsurface  thermocouples  were  inserted, 
o  The  surface  thermocouple  was  embedded  in  a  globule  of  high  conduct¬ 
ivity  silicon  grease. 

o  A  firebrick  hole  cover  was  placed  over  the  core  hole  to  prevent 
heat  from  entering  through  the  open  core  hole, 
o  The  heater  was  lit  and  the  desired  heater  temperature  (T^)  was  set 
by  adjusting  the  propane  pressure  (the  propane  pressure  only 
approximately  related  to  the  heater  temperature). 
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o  When  the  heater  had  reached  steady  state,  the  temperature  measure¬ 
ments  were  begun  to  acquire  the  initial  temperature  profile  of 
the  pavement. 

o  The  heater  was  then  centered  over  the  core  hole  and  temperature 
measurements  were  continued  for  the  predetermined  duration  of 
test. 

Two  types  of  tests  were  conducted  at  Loring  AFB.  One  was  a  "slow"  test 
characterized  by  low  heater  temperatures  (around  300°C)  and  a  long  exposure  time 
(around  25  minutes),  and  the  other  was  a  "quick"  test  characterized  by  high 
heater  temperatures  (around  600°C)  and  a  shorter  exposure  time  (1-3  minutes). 
During  the  "slow"  tests,  the  heater  propane  pressure  was  adjusted  such  that  the 
surface  temperature  of  the  pavement  did  not  rise  above  about  230°C  (close  to 
the  flash  point  of  asphalt).  For  the  "quick"  tests  the  heater  was  removed  when 
the  pavement  surface  reached  about  230°C.  Temperatures  were  recorded  every  20 
seconds  in  the  "slow"  test  and  every  3  seconds  in  the  "quick"  test. 

Figure  I I I -3  shows  the  comparison  of  predicted  and  observed  temperatures 
for  one  of  the  slow  tests  carried  out  at  Loring  Air  Force  Base.  The  predicted 
temperatures  are  very  close  to  the  measured  temperatures  at  points  beneath  the 
pavement  surface  throughout  the  duration  of  the  test.  However,  the  computer 
code  consistently  overestimated  the  measured  temperature  at  the  pavement  sur¬ 
face.  It  is  believed  that  this  was  due  to  conduction  losses  of  the  thermocouple 
leads  at  the  pavement  surface.  Sparrow  (27)  indicates  that  it  is  not  possible 
to  estimate  the  magnitude  of  the  temperature  error  due  to  conduction  losses 
with  any  precision.  However,  the  error  is  in  the  direction  that  was  predicted 
by  the  numerical  model,  i.e.  the  true  surface  temperature  is  higher  than  mea¬ 
sured. 

The  observed  and  predicted  temperatures  at  the  typical  removal  depth  (25 
mm)  are  very  close.  In  this  test  with  a  heater  temperature  of  525°C  the  removal 
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EQUIVALENT 


DIMENSION-ESS  TIME  (Fo*10  ) 

Predicted  Spatial -Temperature  History  for  a  Slow  Test 


temperature  (90  C)  is  not  reached  at  the  typical  25  mm  removal  depth  for  about 
13. 5  minutes. 

Figure  I I 1-4  shows  observed  and  predicted  temperature  histories  for  a  quick 
test.  Again,  the  predicted  surface  temperatures  were  higher  than  observed  be¬ 
cause  of  conduction  losses  as  discussed  above.  The  model  did  an  excellent  job 
of  predicting  the  subsurface  temperatures. 

The  close  agreement  between  observed  and  predicted  temperatures  lends 
credibility  to  the  numerical  model  and  indicates  that  the  model  can  be  used  with 
confidence  to  predicted  in  situ  temperature  due  to  surface  heating  by  fossil 


fuel  heaters. 
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A  PARAMETRIC  STUDY  OF  THE 


HEATING  PROCESS  FOR 
ASPHALTIC  CONCRETE  PAVEMENT 


CHAPTER  IV-A  INTRODUCTION 


Preliminary  Remarks 

Part  IV  presents  a  parametric  study  of  the  heating  process  for  asphaltic 
concrete  pavement  using  the  HEAT1  code.  A  detailed  description  and  a  complete 
listing  of  this  code  appear  in  Part  II  and  Appendix  II-D,  respecti vely. 

A  Typical  Strategy  for  Heating  Asphaltic  Pavement  Prior  to  Removal 
of  a  Surface  Layer. 

The  HEAT1  code  models  the  thermal  aspects  of  a  recycling  process  for  as¬ 
phaltic  concrete  pavement.  Specifically,  the  code  models  that  phase  of  the 
recycling  process  during  which  a  pavement  is  heated  at  its  upper  surface  by  a 
fossil  fuel-fired  heater  prior  to  removal  of  the  upper  part  of  the  pavement  by 
mechanical  means.  The  object  of  the  heating  process  is  to  reduce  the  shearing 
strength  of  the  entire  layer  of  pavement  to  be  removed  by  raising  the  tempera¬ 
ture  of  this  layer  to  above  a  minimum  desirable  temperature  for  removal  (T  , 

removal 

=  363°K)  without  exceeding  a  maximum  allowable  temperature  (=  505°K)  anywhere  in 
the  pavement. 

For  economic  reasons,  the  heating  process  must  be  completed  as  rapidly  as 
possible.  Due  to  the  extremely  low  thermal  diffusivity  of  asphaltic  concrete 
(See  Chapters  I-E  and  I-F)  however,  a  typical  pavement  heater  is  capable  of 
raising  the  surface  temperature  of  the  pavement  to  above  the  maximum  allowable 
value  (505°K)  long  before  the  temperature  at  the  removal  depth  has  reached  the 
minimum  desirable  temperature  for  removal  (363°K).  Thus,  some  contractors  have 
adopted  a  strategy  for  heating  asphaltic  pavement  that  may  be  summarized  in  the 
following  steps: 

(1)  The  pavement  is  heated  until  the  maximum  allowable  temperature  is  reached  at 

the  upper  surface. 
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(2)  Heating  is  then  discontinued  and  the  surface  is  insulated  for  an  arbitrar¬ 
ily  specified  period  of  time  during  which  the  surface  temperature  deceases 
due  to  the  diffusion  of  thermal  energy  into  the  interior  of  the  pavement  by 
conduction. 

(3)  The  insulation  is  removed  from  the  surface  and  the  heating  is  begun  again. 

Ideally,  this  heating  insulation  cycle  is  repeated  until  the  temperature  reaches 
the  desired  value  (363°K)  at  the  depth  which  is  to  be  removed.  From  a  practical 
standpoint,  however,  a  heating  process  that  requires  more  than  two  heating- insu¬ 
lation  cycles  would  be  highly  undesirable  for  use  in  the  field  by  a  contractor. 

CHAPTER  IV-B  DESCRIPTION  OF  THE  PARAMETRIC  STUDY 

Preliminary  Remarks 

HEAT1  uses  a  finite-difference  technique  to  perform  numerical  solutions  of 
the  analytical  model  which  governs  the  heat  transfer  process  in  an  asphaltic  con¬ 
crete  pavement.  The  analytical  model  has  been  recast  in  nondimensional  form, 
thereby  reducing  the  number  of  parameters  from  fourteen  to  nine.  Thus,  the  para¬ 
metric  study  of  the  heating  process  consists  of: 

(1)  defining  a  "nominal  heating  case";  that  is,  selecting  a  set  of  nominal 
values  for  the  nine  dimensionless  parameters  in  the  model; 

(2)  determining  the  solution  of  the  problem  for  the  nominal  heating  case; 

(3)  systematically  varying  the  values  of  the  parameters  in  the  model  to  deter¬ 
mine  their  effects  on  the  heating  process  (specifically,  the  effects  of 
these  variations  on  the  required  heating  time); 

(4)  providing  physical  interpretations  of  the  numerical  results; 
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(5)  and  identifying  (whenever  possible)  combinations  of  the  dimensionless  para¬ 
meters  that  result  in  relatively  short  heating  times. 


Note:  This  is  a  parametric  study,  not  an  optimization  study.  No  attempt 
has  been  made  to  identify  an  optimum  set  of  the  dimensionless  parameters 
which  would  yield  an  absolute  minimum  heating  time. 

i 


Discussion  of  the  Dimensionless  Parameters 

The  fourteen  parameters  in  the  dimensional  analytical  model  of  the  heating 
process  are  defined  below: 


hc  =  convective  heat  transfer  coefficient 

k  =  thermal  conductivity  of  the  pavement 

L  =  overall  thickness  of  the  pavement 
* 

t  TK1{;|||  =  duration  of  the  insulation  period  for  the  upper  surface 
of  the  pavement 

* 

T  H  =  heater  surface  temperature 
* 

T  =  maximum  allowable  temperature  in  the  pavement 
* 

T  p  =  an  appropriate  reference  temperature  for  the  convective  heat 
transfer 

* 

T  =  temperature  to  which  the  pavement  must  be  heated  at  the 

* 

depth  to  which  it  is  to  be  removed  (x  REM) 

* 

T  2qjl  =  temperature  at  the  bottom  surface  of  the  pavement 
* 

x  =  depth  of  pavement  which  is  to  be  removed 
a  s  thermal  diffusivity  of  the  pavement 
=  emissivity  of  the  heater  surface 
=  emissivity  of  the  upper  surface  of  the  pavement 
o  =  Stefan-Boltzman  constant 
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The  minimum,  maximum,  and  nominal  values  of  the  fourteen  dimensional  para 


meters  are  shown  in  Table  IV-1.  The  data  in  Table  I V- 1  show  that  the  values  of 
*  *  * 

the  parmeters  T  T  p^,  x  p^,  and  a  remain  constant  throughout  the  para¬ 

metric  study  and  that  the  value  of  the  reference  temperature  for  convective  heat 

*  * 
transfer  (T  p)  is  assumed  to  be  equal  to  the  heater  surface  temperature  (T  ^). 

Also,  note  that  the  nominal  heating  case  assumes  that  the  heat  transfer  at  the 


upper  surface  of  the  pavement  is  by  radiation  alone  (no  convective  heat  trans¬ 


fer).  Finally,  it  should  be  mentioned  at  this  point  that  the  entire  parametric 


study  (with  two  exceptions  to  be  discussed  later)  was  performed  under  the 

* 

assumption  of  a  uniform  initial  temperature  profile  in  the  pavement  (T 
* 


When  the  analytical  model  is  recast  in  nondimensional  form,  the  fourteen 


dimensional  parameters  defined  above  are  combined  to  give  the  nine  dimensionless 


parameters  shown  below: 


A  n 


a  L  T 


REM 


k(I  +  i  -!) 

£H  £S 


=  radiation  -  conduction  parameter 
=  (thermal  conductance  for  radiation  heat  transfer) 
(thermal  conductance  for  conduction  heat  transfer) 


Bi  =  c 

-  =  Biot  number  = 

k 


(thermal  conductance  for 
convection  heat  transfer) 
(thermal  conductance  for 
conduction  heat  transfer) 


* 

^INSUL  =  t  INSUL  =  dimensionless  duration  of  the  insulation  period 
_  for  the  upper  surface  of  the  pavement 

(L  /or) 


* 

^  H  =  dimensionless  heater  temperature 
- 55 - 

T  REM 


TMAX 


* 


T  MAX 
~H — 

T  REM 


1.391)  =  dimensionless  maximum  allowable  temperature 
in  the  pavement 
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DIMENSIONAL 

MINIMUM 

MAXIMUM 

NOMINAL 

PARAMETER 

VALUE 

VALUE 

VALUE 

hc,  C  — ) 
m2-°k 

0 

300 

0 

k,  (  — - — .  ) 

M-°K 

0.9 

1.65 

1.2 

L,  (M) 

0.254 

0.381 

0.3048 

t  INSUL’  (Sec) 

60 

300 

300 

t*h,  (°K) 

573.16 

1253.16 

673.16 

t*max’  (°k> 

505 

505 

505 

t*r,  (°K) 

573.16 

1253.16 

673.16 

T  REM’  (°K) 

363.16 

363.16 

363.16 

T  SOIL’  (°K) 

273.16 

310.94 

293.16 

x  REM,  (M) 

0.0254 

0.0254 

0.0254 

or.  (  M2  ) 

3.7  x  10“7 

1  x  10-6 

6.6  x  10*7 

Sec 

<-) 

0.3 

0.58 

0.44 

es»  C-) 

0.75 

0.95 

0.9 

o,  (  ) 

5.6697  x  10'8 

5.6697  x  10‘8 

5.6697  x  10‘8 

m2-°k4 

TABLE  IV-1.  Summary  of  the  Minimum,  Maximum,  and  Nominal  Values  of  the  Fourteen 
Parameters  in  vhe  Dimensional  Model  of  the  Heating  Process. 


=  dimensionless  reference  temperature  for  convective 
heat  transfer 


TREM  h  T  REM 


REM 


=  dimensionless  temperature  to  which  the  pavement  must 
be  heated  at  the  depth  to  which  it  is  to  be  removed 


(X 


REM 


) 


t 


TS0IL  e 


* 


T 

S 

T 


SOIL  =  dimensionless  temperature  at  the  bottom  of  the 
pavement  (the  initial  dimensionless  temperature 
REM  throughout  the  pavement) 


* 

v  v 

REM  5  REM  =  dimensionless  depth  of  the  pavement  which  is  to  be 
L  removed. 


The  minimum,  maximum,  and  nominal  values  of  the  nine  dimensionless  para¬ 
meters  are  given  in  Table  IV-2.  These  values  were  computed  using  the  data  in 
Table  IV-1  and  are  the  values  used  to  perform  the  parametric  study. 

Selection  of  the  Dimensionless  Spatial  and  Time  Increments 

The  finite-difference  solution  of  the  heating  process  performed  by  the  HEAT1 
code  requires  the  selection  of  dimensionless  spatial  and  time  increments.  On 
the  basis  of  a  number  of  trial  runs,  the  following  values  were  chosen  for  these 
quantities: 

(dimensionless  spatial  increment)  =  Ax  =  .01  (101  nodal  points) 

(dimensionless  time  increment)  =  At  =  10  ® 

The  entire  parametric  study  was  performed  using  these  values. 
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DIMENSIONLESS 

MINIMUM 

MAXIMUM 

NOMINAL 

PARAMETER 

VALUE 

VALUE 

VALUE 

A 

0.114 

0.647 

0.289 

Bi 

0 

127 

0 

-4 

-3 

-3 

tINSUL 

1.529  x  10 

2.131  x  10 

2.131  x  10  ^ 

th 

1.578 

3.451 

1.854 

tmax 

1.391 

1.391 

1.391 

tr 

1.578 

3.451 

1.854 

trem 

1.0 

1.0 

1.0 

tsoil 

0.752 

0.856 

0.807 

XREM 

0.0667 

0.1 

0.0833 

TABLE  IV-2.  Minimum,  Maximum,  and  Nominal  Values  of  the  Nine 
Dimensionless  Parameters  used  to  Perform  the 
Parametric  Study. 
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CHAPTER  IV-C  RESULTS  OF  THE  PARAMETRIC  STUDY 


The  Nominal  Heating  Case 

The  nominal  heating  case  is  defined  by  the  nominal  values  of  the  nine  dimen¬ 
sionless  parameters  listed  in  Table  IV-2.  Representative  temperature  histories 
and  distributions  are  presented  here  in  graphical  form.  Figure  IV-1  shows  the 
surface  and  removal  depth  temperature  histories.  The  monotonical ly  rising  sur¬ 
face  temperature  clearly  shows  that  no  insulation  periods  are  required  for  the 
nominal  case.  The  heater  temperature  is  low  enough  that  thermal  energy  can 
diffuse  into  the  pavement  by  conduction  fast  enough  to  raise  the  temperature  at 
the  removal  depth  to  its  required  value  (=1.0)  before  the  surface  temperature 

exceeds  the  maximum  allowable  value  (=1.391)  in  the  pavement.  The  Fourier  num- 

-3 

ber  at  the  end  of  the  heating  period  is  (13.4  x  10  ),  which  corresponds  to  a 
heating  time  of  31.5  (min).  From  a  practical  point  of  view,  this  is  an  unac¬ 
ceptably  long  heating  period. 

Figure  IV-2  shows  the  dimensionless  temperature  profiles  in  the  pavement  at 
different  dimensionless  times  during  the  heating  period.  These  profiles  dra¬ 
matically  illustrate  the  exceptionally  low  value  of  the  thermal  diffusivity  of 
asphalt.  At  the  end  of  the  heating  period  (31.5  min),  thermal  energy  from  the 
surface  heating  has  diffused  through  only  about  40%  of  the  total  thickness  of 
the  pavement. 


The  Effect  Of  Initial  Pavement  Temperature  On  Heating  Time 

Figure  IV-3  shows  the  effect  of  the  initial  pavement  temperature  (T$qj|_)  on 
the  heating  time  for  an  otherwise  nominal  case.  The  minimum  value  of  T^^ 

-3 

(0.752)  yields  a  dimensionless  heating  time  of  (17.89  x  10  ),  or  42.05  min.; 

while  the  maximum  value  of  T^g^  (0.856)  requires  a  heating  time  of  (9.822  x 
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Figure  IV-2.  Dimensionless  Temperature  Profiles  In  The 
Pavement  For  The  Nominal  Heating  Case 
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-3 

10  ),  or  23.08  min.  Comparing  these  times  with  the  nominal  heating  time  of 
-3 

(13.4  x  10  ),  or  31.5  min.,  shows  as  expected,  that  initial  temperatures  lower 
or  higher  than  the  nominal  initial  value  (0.807)  significantly  increase  or 
decrease,  respectively,  the  required  heating  time.  Note,  however,  that  even  the 
minimum  heating  time  of  23.08  min.  is  unacceptable  large. 

All  of  the  above  calculations  were  made  by  assuming  a  uniform  initial  tem¬ 
perature  distribution  in  the  pavement.  It  is  well-known,  however,  that  the 
initial  temperature  distribution  in  an  actual  pavement  is  nonuniform,  although 
the  general  functional  form  of  this  distribution  is  unknown.  Thus,  in  order  to 
assess  the  validity  of  performing  the  parametric  study  with  a  uniform  initial 
temperature  distribution,  two  heating  cases  were  computed  using  nonuniform  (1  in- 
near)  initial  temperature  profiles.  These  profiles,  which  are  specified  in  the 
HEAT1  code,  are  based  on  soil  and  surface  temperatures  representative  of  typical 
"hot"  and  "cold"  days.  The  results  of  these  two  cases  are  also  displayed  in 
Figure  IV-3.  The  plotting  of  these  two  points  is  based  on  the  assumption  that 
the  appropriate  values  of  T<-q^  are  given  by  the  average  of  the  soil  and  surface 
temperatures  for  each  profile.  The  heating  time  for  the  linear  initial  tempera¬ 
ture  profile  for  a  hot  day  is  seen  to  be  less  than  the  heating  time  for  a  uni¬ 
form  initial  temperature  distribution  with  the  same  value  of  T<-qjL  (0.8447). 
This  arises  from  the  fact  that  the  initial  temperatures  in  the  layer  to  be  re¬ 
moved  are  higher  than  TSQIL  =  0.8447  for  the  hot  day  case  with  a  linear  initial 
temperature  profile.  Analogous  reasoning  can  be  used  to  explain  why  the  data 
point  for  the  cold  day  heating  case  lies  above  the  (extrapolated)  heating  time 
versus  initial  temperature  curve.  These  calculations  clearly  indicate  that  the 
heating  times  computed  by  assuming  uniform  initial  temperature  profiles  are 
highly  representative  of  the  heating  times  that  would  result  from  the  (more  com¬ 
plicated)  computations  based  on  more  realistic  nonuniform  initial  temperature 
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profiles.  Thus,  all  subsequent  calculations  in  the  parametric  study  are  based 
on  the  assumption  of  uniform  initial  temperature  profiles. 


The  Effects  Of  Heater  Temoerature  and  the  Duration  Of  the  Insulation  Period  on 


Heating  Time 

The  heating  time  is  shown  as  a  function  of  the  heater  temperature  (T^)  in 

Figure  IV-4.  The  parameter  in  Figure  IV-4  is  the  duration  of  the  insulation 

-3 

period,  j^sUL *  assumes  its  maximum  value  (2.131  x  10  ) ,  an  intermediate 

-3  -3 

value  (1.2  x  10  ),  and  its  minimum  value  (0.1529  x  10  ).  The  other  seven 

dimensionless  parameters  retain  their  nominal  values  during  this  portion  of  the 

parametric  study. 

Figure  IV-4  clearly  shows  that,  for  any  given  value  of  the  heating 

time  -  heater  temperature  curve  is  a  U-shaped  curve  which  is  concave  upward. 
Thus,  there  are  two  values  of  the  heater  temperature  which  correspond  to  any 
given  value  of  the  heating  time.  The  entire  curve  for  a  given  value  of  tjN5UL 
is  composed  of  two  distinct  branches:  a  "left-hand  branch"  for  T^  <  2.257  and 
a  "right-hand  branch"  for  T^  >  2.257.  As  shown  in  Figure  IV-4,  the  left-hand 
branch  is  independent  of  the  value  of  tj^suL'  This  branch,  which  represents 
the  locus  of  all  those  operating  points  for  which  no  insulation  periods  are 
required,  includes  the  nominal  operating  point.  Thus,  for  all  of  the  operating 
points  on  the  left-hand  branch  of  the  curve  the  thermal  processes  in  the  pave¬ 
ment  are  qualitatively  the  same  as  those  for  the  nominal  point,  which  are  por¬ 
trayed  in  Figures  IV-1  and  IV-2.  For  all  such  cases  the  temperature  at  the 
removal  depth  reaches  the  desired  value  for  removal,  T^^,  before  the  surface 
temperature  reaches  the  maximum  allowable  value  in  the  pavement,  There¬ 

fore,  when  operating  on  the  left-hand  branch  of  the  curve,  the  heating  time  can 
always  be  reduced  by  increasing  the  heater  temperature.  Such  reductions  are 
possible,  of  course,  because  increases  in  the  heater  temperature  result  in  allow 
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able  increases  in  the  surface  temperature.  These  increases  in  the  surface  tem¬ 
perature  steepen  the  time-averaged  temperature  gradient  which  serves  as  the 
potential  for  the  transfer  of  thermal  energy  by  conduction  in  the  pavement. 


This  behavior  continues  until  reaches  a  value  of  2.257,  which  corresponds  to 
.  .  .  .  -3 

the  minimum  heating  time  of  5.947  x  10  for  which  no  insulation  periods  are 
required.  At  this  point,  the  surface  temperature  reaches  TM^  at  exactly  the 
same  instant  as  the  temperature  at  the  removal  depth  reaches  Tp^.  Any  further 
increases  in  cause  the  surface  temperature  to  reach  before  the  tempera¬ 

ture  at  the  removal  depth  reaches  T^^,  thereby  necessitating  the  use  of  at 
least  one  insulation  period.  Thus,  the  right-hand  branch  of  the  heating  time 
versus  heater  temperature  curve  corresponds  to  the  locus  of  all  those  operating 
points  for  which  one  or  more  insulation  periods  are  required.  The  numbers  of 
heating-insulation  cycles  for  a  selected  group  of  such  operating  points  are 
displayed  in  Figure  IV-4. 

At  first  glance,  the  right-hand  branch  of  the  curve  exhibits  a  rather  un¬ 
expected  behavior;  that  is,  increases  in  heater  temperature  result  in  correspond¬ 
ing  increases  in  the  required  heating  time.  Qualitatively  this  behavior  is 
independent  of  the  value  of  tj^uL’  however,  as  may  be  seen  in  Figure  IV-4,  the 


magnitude  of  these  increases  in  heating  time  is  a  strong  function  of  t 


INSUL  ’ 


Thus  two  important  features  of  the  right-hand  branch  of  the  curve  must  be  ex¬ 
plained:  (1)  the  qualitative  effects  of  ^ INSUL  on  the  Seating  time  ^or  a  9iven 
heater  temperature,  and  (2)  the  relative  independence  of  the  heating  time  on  the 


heater  temperature  for  small  values  of  t 


INSUL’ 


As  shown  in  Figure  IV-4,  when  t 


INSUL 


assumes  its  largest  possible  value 


(2.131  x  10  J),  the  heating  time  varies  over  a  wide  range  of  values  as  a  func¬ 


tion  of  the  heater  temperature.  This  trend  for  a  large  value  of  t 


INSUL 


easily  explained  with  the  aid  of  Figure  IV- 5 ,  which  is  a  plot  of  the  surface 
and  removal  depth  temperature  histories  for  a  heater  temperature  of  2.75  and 
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tjNsui  ~  2.131  x  10  At  heater  temperatures  greater  than  2.2b 7,  the  surface 
temperature  rises  quickly  and  reaches  before  the  temperature  at  the  removal 
depth  reaches  Tp^.  As  may  be  seen  in  Figure  IV-5,  this  initiates  an  insulation 
period  of  a  duration  which  is  significantly  longer  than  that  of  a  typical  heat¬ 
ing  period.  During  each  of  these  "long"  insulation  periods,  thermal  energy  dif¬ 
fuses  into  the  interior  regions  of  the  pavement  by  conduction;  however,  the  flow 
of  additional  thermal  energy  into  the  pave*  'nt.  at  its  upper  surface  is  delayed 
until  the  end  of  the  insulation  period.  Thus,  the  net  effect  of  employing  a 
long  insulation  period  is  to  markedly  flatten  out  the  temperature  gradient  in 
the  pavement.  This  reduces  the  time-averaged  rate  of  conduction  heat  transfer 
into  the  interior  of  the  pavement  and  thereby  increases  the  required  heating 
time. 


In  contrast  to  this  behavior  for  a  large  value  of  tjNSUL.  Figure  IV-4  re¬ 
veals  that  the  behavior  of  the  right-hand  branch  is  qual i tati vely  different  for 
small  values  of  ^ I NSUL ’  fhlS  d^frerence  behavior  may  be  best  explained  in 
terms  of  an  optimal  (minimum  time)  heating  process.  In  this  idealized  optimal 
heating  process  the  initial  heater  temperature  would  be  so  large  that  the  sur¬ 
face  temperature  would  rise  to  T^^  in  a  very  short  period  of  time;  however, 
this  would  not  initiate  an  insulation  period.  Instead,  the  heater  temperature 
would  be  reduced  in  a  manner  such  that  the  surface  temperature  was  maintained  at 
a  constant  value  of  until  the  temperature  at  the  removal  depth  had  risen  to 
TreM'  process  would  result  in  the  maximum  possible  time-averaged  tempera¬ 
ture  gradient  in  the  pavement  and,  therefore,  in  the  minimum  possible  heating 
time.  (The  absolute  minimum  heating  time  would  be  obtained  in  the  limiting 
case  of  the  optimal  heating  process  which  corresponds  to  a  step  change  in  the 


surface  temperature  from 
required  heating  time  in  this 


to  Under  otherwise  nominal  conditions,  the 

limit  is  (3.66  x  10  ^);  or  approximately  8.6  min¬ 


utes).  Implementation  of  the  optimal  heating  process  obviously  requires  the  use 


of  a  variable  temperature  heater;  however,  this  is  forbidden  by  the  scope  of  the 
present  study.  It  is  possible,  though,  to  operate  a  constant  temperature  heater 
such  that  the  deviation  of  the  actual  surface  temperature  history  from  the 
optimal  surface  temperature  history  is  arbitral' ly  small.  This  can  be  accom¬ 
plished  by  decreasing  the  duration  of  the  insulation  period,  tj^<-^.  For  3  very 
small  value  of  tj^^  there  is  only  a  slight  decrease  in  the  surface  temperature 
from  during  an  insulation  period,  after  which  the  heater  quickly  raises  this 

temperature  back  to 

Thus,  in  the  limit  of  tjNSU^  approaching  zero  the  surface  temperature 
oscillates  in  a  very  narrow  band  of  temperatures  with  an  upper  bound  of  T^^. 
The  width  of  this  band  can  always  be  made  less  than  any  desired  value  by  choos¬ 
ing  a  small  enough  value  of  t^yy  Since  the  duration  of  a  typical  heating 
period  also  approaches  zero  in  the  limit  of  tJNSUL  approaching  zero,  the  surface 
temperature  history  very  closely  approximates  the  optimal  surface  temperature 
history  which  accounts  for  the  low  heating  times  for  small  values  of  tjNSUL. 
Furthermore,  the  difference  between  the  surface  temperature  at  the  end  of  an 
insulation  period  and  T^^  approaches  zero  as  tj^yy  approaches  zero.  Since 
this  temperature  difference  is  small  and  the  value  of  T^  is  relatively  large, 
the  duration  of  the  ensuing  heating  period  will  be  essentially  independent  of 
the  value  of  T^.  Thus,  since  ^j^gyy  1s  3  fixed  value,  the  required  heating  time 
will  be  a  weak  function  of  the  heater  temperature. 

The  behavior  described  above  for  small  values  of  tjN^yy  is  portrayed  in 

Figures  IV-4  and  IV-6.  In  Figure  IV-4  the  right-hand  branch  of  the  heating  time 

-3 

versus  heater  temperature  curve  for  the  minimum  value  of  j NSUL  (0-1529  x  10  ) 

exhibits  the  expected  relatively  low  values  of  the  heating  time  and  the  weak 
dependence  of  the  heating  time  on  heater  temperature.  Figure  IV-6  is  a  plot  of 


the  surface  and  removal  depth  temperature  histories  for  a  heater  temperature  of 

-3 

2.75  and  the  shortest  possible  insulation  period  (0.1529  x  10  ).  [Recall  that 
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Figure  IV-5  is  a  plot  of  these  same  temperature  histories  for  Tu  =  2.75  and  the 
maximum  value  of  tj^^y  (2.131  x  10  ^)].  Figures  IV-4,  IV-5,  and  I V- 6  illus¬ 
trate  the  effect  of  tj^yy  on  the  heating  time,  which  is  reduced  from  (9.845  x 
-3  -3 

10  )  to  (4.688  x  10  )  as  ^ I NSU L  ^ecreases  ^rom  maximum  value  of  (2.131  x 

~3  -  3 

10  )  to  its  minimum  value  (0.1529  x  10  )  at  a  heater  temperature  of  2.75. 

Figures  IV-5  and  IV-6  also  portray  the  approach  of  the  surface  temperature  his¬ 
tory  to  the  optimum  surface  temperature  history  in  the  limit  of  tjNSyL  approach¬ 
ing  zero  as  described  above. 

The  forgoing  discussions  might  lead  to  the  erroneous  conclusion  that  the 

best  heating  process  for  an  asphaltic  concrete  pavement  is  one  which  employes  a 

heater  temperature  greater  than  2.257  and  the  smallest  possible  value  of  t^^y. 

In  the  present  study  this  would  reduce  the  heating  time  from  the  minimum  value 
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required  for  no  insulation  periods  (13.4  x  10  )  to  the  minimum  value  of  (4.688 

-3 

x  10  )  corresponding  to  a  heater  temperature  of  2.75  at  the  smallest  possible 

-3 

value  of  tj^SUy  (0.1529  x  10  ),  a  reduction  of  65  percent.  As  shown  in  Figure 
IV-6,  however,  the  number  of  heating-insulation  cycles  [14]  required  to  achieve 
this  reduction  is  unacceptably  large  for  practical  implementation. 

Based  on  the  above  discussion,  the  following  conclusions  may  now  be  made: 

(1)  The  minimum  possible  heating  time  may  be  achieved  with  heater  temperatures 
greater  than  2.257  and  with  the  minimum  value  of  ^ jj^suy ' 

(2)  This  is  not  a  practical  mode  of  operation,  however,  since  the  number  of 
heating-insulation  cycles  is  prohibitively  large. 

(3)  The  use  of  the  "heating-insulation  cycle  strategy"  by  paving  contractors 
should  be  discontinued  since  it  is  not  a  practical  technique. 

(4)  Paving  contractors  who  wish  to  employ  surface  heating  techniques  should 
operate  their  heaters  at  the  maximum  possible  temperature  (2.257  in  this 


study)  below  which  no  heating-insulation  cycles  are  required. 

Note:  This  recommendation  is  based  on  the  assumption  of  no  convective  heat 
transfer.  As  will  be  shown,  a  more  economical  (lower  than  2.257)  heater 
temperature  may  be  employed  if  there  is  some  convective  heating  at  the  upper 
surface  of  the  pavement. 

The  Effect  of  the  Biot  Number  on  the  Heating  Time 

The  Biot  number,  which  is  defi.ied  by  the  expression  (hcL/k),  may  be  inter¬ 
preted  physically  as  the  ratio  of  the  thermal  conductance  for  convective  heat 
transfer  to  the  thermal  conductance  for  conduction  heat  transfer.  For  the 
nominal  heating  case  all  of  the  heat  transfer  at  the  upper  surface  of  the  pave¬ 
ment  is  assumed  to  take  place  by  radiation.  Thus,  the  Biot  number  for  this  case 
is  identically  equal  to  zero.  As  the  Biot  number  increases  from  zero,  the  net 
rate  of  heat  transfer  to  the  pavement  (radiation  plus  convection)  increases  due 
to  the  increase  in  convection  heat  transfer  with  no  corresponding  increase  in 
the  heater  temperature.  Therefore,  qualitatively,  an  increase  in  the  Biot  num¬ 
ber  should  have  the  same  effect  on  the  required  heating  time  as  an  increase  in 
the  heater  temperature.  This  behavior  is,  indeed,  displayed  in  Figure  IV-7, 
which  shows  the  required  heating  time  as  a  function  of  the  Biot  number  for  an 
otherwise  nominal  heating  case.  The  characteristic  U-shaped  curve  may  be  inter¬ 
preted  in  the  same  manner  as  the  heating  time  versus  heater  temperature  curve 
(Figure  IV-4).  Of  special  interest  here  is  the  left-hand  branch  of  the  curve 
which  corresponds  to  the  locus  of  all  those  operating  points  for  which  no  heat¬ 
ing-insulation  cycles  are  required.  When  operating  on  this  branch,  increases  in 
the  Biot  number  result  in  decreases  in  the  required  heating  time.  This  behavior 

persists  up  to  a  Biot  number  of  10,  which  corresponds  to  a  dimensionless  heating 
-3 

time  of  (4.694  x  10  ),  or  11.75  minutes  under  nominal  conditions. 
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Since  heater  fuel  consumption  is  directly  proportional  to  some  function  of 

the  heater  temperature,  it  is  always  desirable  to  operate  the  heater  at  the 

lowest  possible  temperature.  Thus,  an  important  result  may  now  be  found  by  not- 

-3 

ing  that  the  minimum  heating  time  of  (4.964  x  10  )  for  an  off-nominal  Biot  num¬ 
ber  of  10  was  computed  for  a  nominal  heater  temperature  of  1.854.  This  time 

-3 

corresponds  to  a  16.6%  reduction  from  the  minimum  heating  time  of  (5.947  x  10  ) 

computed  for  a  nominal  value  (zero)  of  the  Biot  number  and  an  off-nominal  heater 
temperature  of  2.257. 

Therefore,  we  may  conclude  that: 

(1)  The  use  of  convective  heating  is  desirable. 

(2)  Operating  points  on  the  right-hand  branch  of  the  heating 
time  versus  Biot  number  curve  are  always  impractical. 

The  first  conclusion  is,  of  course,  based  on  the  fact  that  the  minimum  heating 

-3 

time  may  be  reduced  to  (4.964  x  10  )  by  employing  convective  heating  such  that 

the  Biot  number  is  10.  (Additional  information  on  the  desirability  of  convective, 
as  opposed  to  radative,  heating  will  be  included  in  the  following  section).  The 
second  conclusion  follows  since  operating  points  on  the  right-hand  branch  of  the 
curve  correspond  to  longer  than  minimum  heating  times  and  require  prohibitively 
large  number  of  heating-insulation  cycles  as  shown  in  Figure  IV-7. 

The  Effect  of  the  Radiation-Conduction  Parameter  ("A")  on  the  Heating  Time 

Figure  IV-8  displays  the  effect  of  the  radiation-conduction  parameter, 
which  is  defined  by  the  expression: 

i  t*3 
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on  the  heating  time  for  an  otherwise  nominal  heating  case.  This  parameter  may 
be  interpreted  physically  as  the  ratio  of  the  thermal  conductance  for  radiation 
heat  transfer  to  the  thermal  conductance  for  conduction  heat  transfer. 

Over  the  range  of  values  of  A  considered  in  the  present  study  (0.114  - 
0.647),  the  heating  time  is  seen  to  be  a  monotonical ly  decreasing  function  of  A. 
Since  the  Biot  number  is  equal  to  zero  for  the  nominal  heating  case  (no  convec¬ 
tive  heat  transfer)  and  the  parameter  A  is  always  less  than  unity,  the  behavior 
exhibited  in  Figure  IV-8  clearly  indicates  that  the  heating  process  with  no  con¬ 
vection  is  "radiation  controlled".  Thus,  the  overall  rate  of  transfer  of  ther¬ 
mal  energy  from  the  heater  to  the  inner  regions  of  the  pavement  is  limited  in 
this  case  by  the  radiative  heat  transfer  process  at  the  upper  surface  of  the 
pavement,  rather  than  by  the  diffusion  of  thermal  energy  by  conduction  in  the 
pavement.  Therefore,  when  A  increases  (meaning  that  the  thermal  conductance  for 
radiation  increases  relative  to  the  thermal  conductance  for  conduction)  there 
should  be  a  corresponding  decrease  in  the  heating  time  as  shown  in  Figure  IV-8. 
Given  the  extremely  low  thermal  diffusivity  of  asphaltic  concrete,  this  is  an 
interesting  and,  perhaps,  unexpected  result  which  suggests  that  the  use  of  a 
thermal  radiation  heater  is  a  rather  inefficient  technique  for  this  application. 

A  deeper  understanding  of  the  above  result  and  the  result  of  the  previous 
section  can  be  obtained  by  examining  the  relative  importance  of  the  radiative 
and  convective  heat  transfer  process  for  the  otherwise  nominal  case  for  which 
the  Biot  number  is  10.  Keeping  in  mind  the  physical  interpretations  of  the  Biot 
number  (Bi)  and  the  radiation-conduction  parameter  (A),  we  have  for  this  case: 
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thermal  conductance  for  convection  heat  transfer 
thermal  conductance  for  radiation  heat  transfer 
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or,  Bi  =  10  =  34.6; 
A  .289 
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that  is,  the  thermal  conductance  for  convective  heat  transfer  is  approximately 
35  times  greater  than  the  thermal  conductance  for  radiation  heat  transfer.  This 
suggests  that  a  convection  heater,  rather  than  a  thermal  radiation  heater, 
should  be  used  for  the  surface  heating  of  asphaltic  concrete  pavements.  Finally, 
since  the  radiative  heating  is  essentially  negligible  with  respect  to  the  convec¬ 
tive  heating  for  this  case,  the  physical  interpretation  of  the  Biot  number  as 
the  ratio  of  the  thermal  conductance  for  convection  to  the  thermal  conductance 
for  conduction  shows  that,  when  convection  heating  dominates  radiative  heating, 
the  overall  heat  transfer  process  is  "diffusion  limited".  This  means  that  the 
overall  rate  of  transfer  of  thermal  energy  from  the  heater  to  the  inner  regions 
of  the  pavement  is  controlled  by  the  diffusion  of  thermal  energy  by  conduction 
in  the  pavement. 

The  Effect  of  Removal  Depth  on  Heating  Time 

The  effect  of  removal  depth  on  heating  is  displayed  in  Figure  IV-9  for  an 
otherwise  nominal  case.  As  expected,  this  plot  simply  shows  that  increasing  the 
removal  depth  increases  the  required  heating  time. 
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PART  V 

SUMMARY  OF  RESULTS 
AND  CONCLUSIONS 


CHAPTER  V- A  SUMMARY  OF  PART  I-THERMAL  PROPERTIES 


As  a  result  of  the  research  described  Part  I,  the  following  results  and 

conclusions  can  be  stated: 

1.  In  analyzing  heat  transfer  in  asphaltic  concrete  pavement  having  limestone 
or  granite-gneiss  aggregate,  this  study  indicates  that  an  average  value  of 
specific  heat  and  an  average  value  of  diffusivity  can  be  used  which  are 
independent  of  asphalt  content  but  depend  on  the  aggregate  gradation  and 
type.  The  conductivity  used  in  analysis  must  take  into  account  the  asphalt 
content  in  addition  to  the  aggregate  type  and  gradation. 

2.  Both  aggregate  type  and  gradation  have  a  greater  effect  on  the  thermal 
properties  of  asphaltic  concrete  pavement  than  does  the  percentage  of 
asphalt  present  in  the  mix. 

3.  There  are  two  opposing  mechanisms  that  influence  the  thermal  conductivity 
of  an  asphaltic  concrete  mix  as  the  asphalt  content  is  increased.  On  one 
hand  the  replacement  of  air  in  the  voids  by  asphalt  tends  to  increase  the 
conductivity  of  the  mix  because  the  conductivity  of  asphalt  is  much  higher 
than  that  of  air;  on  the  other  hand  the  additional  asphalt  increases  the 
thickness  of  the  coating  around  aggregates  which  tends  to  decrease  the  con¬ 
ductivity  of  the  mix  because  the  conductivity  of  asphalt  is  much  less  than 
that  of  the  aggregate.  The  dominant  mechanism  depends  on  the  asphalt  con¬ 
tent  and  the  mix  properties  such  as  the  largest  particle  sizes  present,  the 
gradation,  and  the  specific  surface. 

4.  For  a  given  aggregate  type,  as  the  percentage  of  air  voids  in  an  asphaltic 
concrete  is  increased,  the  thermal  conductivity  of  the  mix  will  decrease 


106 


because  the  thermal  conductivity  of  air  is  less  than  that  of  asphalt  or 
aggregate.  Also,  entrapped  air  in  the  aggregate  particles  tends  to  de¬ 
crease  the  thermal  conductivity  of  the  mix. 

5.  The  specific  heat  capacity  of  the  six  laboratory  mixes  tested  did  not  vary 
systematically  as  the  asphalt  content  was  increased.  The  specific  heat  is 
more  dependent  on  the  aggregate  gradation  and  mineralogy  than  the  asphalt 
content.  Based  on  these  data,  it  seems  that  an  average  value  of  c  for  each 
mix,  independent  of  asphalt  content,  could  be  used. 

6.  The  thermal  diffusivity  of  the  open  graded  surface  course  mixes  (limestone 
and  granite-gneiss  aggregates)  is  essentially  independent  of  asphalt  con¬ 
tent;  the  diffusivity  of  the  limestone  base  course  mix  and  the  limestone 
and  granite-gneiss  dense  graded  surface  course  mixes  varied  about  12%  over 
the  range  of  asphalt  contents  used  in  the  mixes.  The  diffusivity  of  the 
dense  graded  surface  course  mix  having  lightweight  aggregate  is  much  less 
than  that  of  the  dense  graded  limestone  mix  and  is  essentially  independent 
of  asphalt  content. 

7.  The  thermal  and  physical  properties  analyses  of  the  field  specimens  from 
Myrtle  Beach  and  Loring  Air  Force  Bases  indicate  no  relationship  between 
the  properties  tested  (p,  k,  c,  and  a)  and  asphalt  content.  There  is  a 
relationship  between  the  thermal  properties  (k,  c,  and  a)  and  the  aggregate 
type  present  in  the  samples.  The  average  (of  21  samples)  thermal  conduc¬ 
tivity  of  the  Loring  samples  is  slightly  higher  that  that  of  the  Myrtle 
Beach  samples  due  to  the  lower  conductivity  of  the  quartz  sand  fine  aggre¬ 
gate  present  in  the  Myrtle  Beach  pavement.  The  average  specific  heat  of 
the  Myrtle  Beach  samples  is  much  less  than  that  of  the  Loring  samples, 
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because  the  Loring  slate-basalt  aggregate  has  higher  specific  heat  than  the 
Myrtle  Beach  granite  and  quartz  sand  aggregate. 

8.  The  thermal  properties  (k,  c,  and  a)  of  asphaltic  pavement  samples  can  be 
determined  with  about  5.8%  or  less  uncertainty  by  using  the  apparatus  and 
procedure  described  in  Chapter  I-C  and  Appendix  I-A. 

CHAPTER  V-B  SUMMARY  OF  PART  II  -  NUMERICAL  MODEL 

As  a  result  of  the  research  performed  in  Part  II,  the  following  conclusions 

and  results  can  be  summarized: 

1.  The  errors  incurred  by  modeling  the  process  of  heat  transfer  through  a 

three  layer  asphaltic  concrete  pavement  as  that  of  the  transfer  of  heat 

through  a  one-layer  (homogeneous)  body  are  negligibly  small. 

2.  The  nondimensional  form  of  the  numerical  model  was  programmed  for  execution 
on  a  digital  computer.  In  the  program,  heating  is  begun  and  continues 
until  either  the  surface  temperature  reaches  a  maximum  allowable  value 
(flash  point  of  asphalt),  or  the  temperature  at  the  prescribed  removal 
depth  attains  a  desired  value,  at  which  time  the  process  is  terminated. 
If  the  temperature  at  the  removal  depth  has  not  reached  the  desired  value 
before  the  maximum  allowable  surface  temperature  is  reached,  the  heating  is 
terminated  and  the  surface  of  the  pavement  in  insulated  for  a  specified 
period  of  time  and  the  cycle  is  repeated  until  the  desired  temperature  at 
the  removal  depth  is  obtained. 

3.  A  "nominal"  heating  case  was  described  in  which  (a)  typical  values  of  the 

parameters  in  the  model  were  assumed;  (b)  the  surface  of  the  pavement  was 
assumed  to  be  heated  by  radiation  alone  (no  convection)  and;  (c)  the  ini¬ 
tial  temperature  distribution  in  the  pavement  was  assumed  to  be  uniform. 
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4. 


The  analysis  showed  that  for  the  nominal  heating  case  no  insulation  periods 
were  required  i.e.,  the  heater  temperature  was  low  enough  that  thermal 
energy  could  diffuse  into  the  pavement  by  conduction  fast  enough  to  raise 
the  temperature  at  the  prescribed  removal  depth  to  its  required  value  be¬ 
fore  che  surface  temperature  exceeded  its  maximum  allowable  value.  For  the 
nominal  case,  the  required  heating  time  was  an  unacceptably  long  31.5  min¬ 
utes. 

5.  Because  of  the  very  low  diffusivity  of  asphaltic  concrete,  for  the  heating 

time  of  31.5  minutes  in  the  nominal  case,  only  the  upper  40%  of  the  total 
thickness  of  the  pavement  experienced  an  increase  in  temperature  due  to  the 
surface  heating. 

CHAPTER  V-C  SUMMARY  OF  PART  III  -  VALIDATION  OF  THE  MODEL 

As  a  result  of  comparing  measured  pavement  temperature  histories  from  field 
tests  carried  out  at  Myrtle  Beach  and  Loring  Air  Force  Bases  and  histories  pre¬ 
dicted  for  these  tests  by  the  thermal  model,  the  following  points  can  be  made: 

1.  The  results  of  the  tests  carried  out  at  Myrtle  Beach  Air  Force  Base  in 

March  of  1983  were  considered  to  be  unreliable  because  the  coastal  winds 

prevailing  at  the  time  disturbed  the  temperature  patterns  in  the  air  be¬ 
tween  the  heater  and  pavement  surface.  Because  of  the  wind  gusts,  unpre¬ 
dictable  convective  cooling  of  the  pavement  surface  occurred  and  surface 
temperatures  as  well  as  heater  temperatures  fluctuated  widely.  Thus  no 
meaningful  analysis  of  these  tests  could  be  made. 

2.  The  wind  problem  was  alleviated  during  the  field  tests  carried  out  at 
Loring  Air  Force  Base  in  June  of  1983  by  erecting  a  corrugated  steel  "wind" 
fence  about  3m  square  and  0.8m  high  around  each  test  site. 
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3.  Analyses  of  the  Loring  test  data  indicated  that  the  analytical  model  con¬ 
sistently  predicted  higher  temperature  at  the  pavement  surface  than  were 
measured  by  thermocouples.  This  was  to  be  expected  due  to  conduction  cool¬ 
ing  of  the  thermocouples  beads.  The  error  introduced  in  the  measurements 
could  not  be  calculated  but  the  error  analysis  is  in  the  right  direction 
i.e.  it  indicates  that  the  measured  surface  temperatures  were  lower  than 
what  actually  existed. 

4.  Comparison  of  the  predicted  and  measured  temperature  histories  for  points 
in  the  pavement  below  the  surface  indicated  that  the  model  did  a  fair  to 
excellent  job  in  predicting  temperature  histories.  Generally,  the  pre¬ 
dicted-measured  temperature  data  matched  more  closely  as  the  depth  from  the 
pavement  surface  increased.  At  the  typical  removal  depth  of  2.5  cm,  the 
computer  model  predicted  the  temperature  time  data  extremely  well. 

5.  Although  the  commercially  available  heater  used  in  the  study  was  purported 
to  be  a  radiant  heater,  test  data  tend  to  support  the  belief  that  convec¬ 
tive  cooling  around  the  edges  of  the  heater  was  an  important  characteristic 
and  should  be  incorporated  in  the  analysis. 

CHAPTER  V-D  SUMMARY  OF  PART  IV  -  PARAMETRIC  STUDY 

In  the  parametric  study  a  nominal  heating  case  was  defined  and  the  solution 
of  the  problem  was  obtained  as  described  previously.  Then  the  values  of  the 
parameter  in  the  model  were  systematically  varied  to  determine  their  effects  on 
the  heating  process.  As  a  result  of  this  part  of  the  study,  the  following  con¬ 
clusions  and  results  can  be  stated: 

1.  The  initial  temperature  of  the  pavement  significantly  affects  the  required 
heating  time.  Also,  it  was  found  that  the  heating  time  computed  by  assum- 
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ing  uniform  initial  pavement  temperature  profiles  was  an  excellent  approx¬ 
imation  of  the  heating  time  resulting  from  a  more  realistic  nonuniform 
initial  temperature  profile. 

2.  The  heating  time  -  heater  temperature  curve  was  found  to  be  a  U-shaped 
curve  which  is  concave  upward.  The  bottom  of  the  curve  defines  the  minimum 
heating  time  corresponding  to  the  heater  temperature  at  which  the  surface 
temperature  reaches  its  maximum  allowable  value  at  the  same  time  the  de¬ 
sired  removal  temperature  is  reached  at  the  prescribed  removal  depth. 

3.  For  heater  temperatures  greater  than  that  corresponding  to  the  bottom  of 
the  U-shaped  curve,  the  surface  temperature  reaches  its  maximum  allowable 
value  before  the  removal  temperature  is  reached  at  the  removal  depth.  Thus, 
insulation  periods  are  required  in  this  (the  right  hand)  part  of  the  curve. 

4.  In  the  left  hand  part  of  the  curve  the  heating  time  can  be  decreased  by 
increasing  the  heater  temperature.  In  this  part  of  the  heating  time  - 
heater  temperature  curve,  the  removal  temperature  is  reached  at  the  re¬ 
moval  depth  before  the  surface  temperature  reaches  its  maximum  allowable 
value  and  no  insulation  periods  are  required. 

5.  The  analysis  indicated  that  heating  times  could  be  decreased  below  the 
value  corresponding  to  "no  insulation"  only  if  many  heating  -  insulation 
cycles  of  very  short  duration  were  used.  Since  this  is  not  a  practical 
alternative  to  use  in  the  field,  an  important  finding  of  this  research  is 
that  "soaking"  or  insulation  periods  should  not  be  used  in  heating  as¬ 
phaltic  pavements.  The  best  strategy  is  to  use  the  maximum  temperature  be¬ 
low  which  no  heating  -  insulation  cycles  are  required. 
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Some  convection  in  addition  to  radiant  heat  transfer  from  the  heater  to 
the  pavement  surface  is  desirable  because  it  decreases  the  minimum  required 
heating  time.  Alternatively,  for  a  given  heating  time,  the  required  heater 
temperature  can  be  reduced  if  some  convection  heat  transfer  occurs.  The 
minimum  heating  time  with  no  insulation  cycles  corresponded  to  a  Biot  num¬ 
ber  of  10. 

It  was  found  that  if  there  were  no  convective  heat  transfer,  the  overall 
rate  of  transfer  of  thermal  energy  from  the  heater  to  the  inner  regions  of 
the  pavement  is  limited  by  the  radiative  heat  transfer  process  at  the  sur¬ 
face  of  the  pavement  rather  than  by  diffusion  of  thermal  energy  by  con¬ 
duction  in  the  pavement.  This  suggests  that  use  of  a  radiative  heater  is  a 
rather  inefficient  technique  for  heating  pavements  even  though  the  diffus- 
ivity  of  asphaltic  concrete  is  small. 

Further  analysis  suggested  that  a  combined  convection-thermal  radiation 
heater,  rather  than  a  thermal  radiation  heater  alone,  should  be  used  for 
the  surface  heating  of  asphaltic  concrete  pavement. 
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APPENDIX  I-A 


TEMPERATURE  MEASURING  SYSTEM 

The  TEMPSENSE  Temperature  Monitoring  System  marketed  by  Interactive  Micro¬ 
ware  (29)  in  conjunction  with  an  Apple  11+  48K  computer  was  used  to  collect 
temperature  and  time  data.  This  Appendix  provides  a  brief  description  of  the 
TEMPSENSE  system  and  of  the  hardware  used  in  collecting  the  data. 

The  hardware  provided  with  the  TEMPSENSE  system  includes  an  ADALAB  data 
acquisition  interface  card,  an  ADA-MUX  8-channel  multiplexer,  and  an  ADA-AMP 
variable-gain  instrumentation  amplifier.  The  voltage  from  each  thermocouple 
goes  first  through  the  multiplexer,  which  is  a  switching  device,  then  to  the 
amplifier  where  the  voltage  is  amplified.  The  amplified  voltage  is  then  routed 
to  the  analog  to  digital  converter  on  the  ADALAB  interface  card.  Interpolation 
tables  for  many  types  of  thermocouples  are  stored  in  the  software. 

Copper-constantan  thermocouples  of  size  30  AWG  were  used  to  determine  the 
temperature  at  the  top  of  the  pavement  sample  (T^.),  while  size  24  AWG  thermo¬ 
couples  were  used  for  the  bottom  of  the  pavement  sample  (Tb).  The  TEMPSENSE 
software  guides  the  user  through  adjustment  of  the  amplifier,  calibration  of  the 
thermocouples  using  an  ice  bath  and  another  known  temperature  bath,  and  selec¬ 
tion  of  the  cold  junction  compensation  method.  A  solid  state  sensor  was  used  to 
measure  the  temperature  of  the  cold  junction  (multiplexer  terminal). 

After  the  amplifier  is  properly  adjusted  and  the  thermocouples  are  cali¬ 
brated,  temperature  measurements  may  be  taken.  The  TEMPSENSE  program  prompts 
the  user  for  the  frequency  of  sampling  (up  to  3  samples  per  second)  and  the 
number  of  measurements  to  be  averaged  for  each  reading.  Then,  the  total  dur¬ 
ation  of  sampling  and  the  name  of  the  disk  file  is  selected.  The  time  of  each 
measurement  is  listed  on  the  screen,  along  with  the  temperature  measured  by  each 
thermocouple.  When  the  storage  buffer  is  full,  the  data  are  automatically 
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stored  on  a  floppy  disk,  and  measurements  continue  until  the  end  of  the  pre¬ 
selected  duration  time.  The  data  are  automatically  stored  in  binary  form,  but 
the  program  allows  the  user  to  change  the  binary  file  to  a  text  file  which  may 
be  more  compatible  to  some  users. 
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APPENDIX  I-B 


UNCERTAINTY  ANALYSIS 


This  section  presents  the  precision  of  each  type  of  measurement  (e.g. 
length,  temperature,  etc.)  and  an  uncertainty  analysis  of  each  of  the  thermal 
properties  (as  suggested  by  Kline  and  McClintock,  Ref.  20)  using  the  apparatus 
and  procedure  described  in  Chapter  I-C.  Also  included  is  an  analysis  of  the 
amount  of  heat  loss  from  the  calorimeter  and  from  the  sample  holder.  The  re¬ 
sults  are  summarized  in  Figure  I - Bl. 

Figure  I-Bl  shows  a  schematic  drawing  of  the  calorimeter  and  the  sample 
holder.  A  steady-state,  one-dimensional  heat  transfer  analysis  was  conducted  to 
determine  the  approximate  amount  of  heat  loss  from  the  calorimeter  and  the  sam¬ 
ple  holder.  The  thermal  properties  of  the  insulating  materials  were  obtained 
from  References  1  and  21.  The  analysis  shows  a  heat  loss  of  about  10.7  watts 
from  the  calorimeter  and  about  2.0  watts  from  the  sample  holder.  For  a  typical 
test  duration  of  about  25  minutes,  the  total  heat  loss  is  about  15,910  J  from 
the  calorimeter  and  about  3,010  0  from  the  sample  holder. 

The  percentage  of  uncertainty  of  each  of  the  thermal  properties  was  deter¬ 
mined  in  the  following  manner: 

1.  An  actual  test  data  set  was  chosen  (because  of  its  similarity  to  sev¬ 
eral  other  test  data)  to  represent  a  typical  test  data  set. 

2.  The  uncertainty  of  each  of  the  primary  experimental  measurements  is 
specificed  (called  w  ,  w.,...,  w  ,  where  the  subscript  refers  to  tne 

3D  n 

measurement  and  w  is  the  uncertainty  of  that  measurement). 

3.  Each  thermal  property  equation  (Equations  7,  8,  and  9  of  Chapter  I-C) 


is  partially  differentiated  with  respect  to  each  of  the  primary  experi- 

dX  dX  dX 

mental  measurements  (  r-  •zr  r-,  where  X  is  the  thermal  proper- 

da  ,  dD  ,  .  .  .  ,  dH 

ty  equation  and  a,  b,...,n  refer  to  the  measurement). 
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Figure  I-Bl.  Calculations  of  Potential  Heat  Loss  from  the 
Calorimeter  and  from  the  Sample  Holder. 


4.  Kline  and  McClintock  (20)  have  shown  that  the  uncertainty  (w^)  of  each 
of  the  thermal  properties  can  be  determined  by  the  following  equation: 


it  3X  \2  a  /  OX  .  .  /  OX  \2l^ 

w- =  k  ai  wa  y  (  hk  w..  y 2  +  (  ^  > 2i 


(I-Bl) 


3b  b 


3n  n 


where  w^  =  uncertainty  of  the  value  of  equation  X 

wn  =  uncertainty  of  the  primary  measurement  n 
3  X 

—  =  partial  of  equation  X  with  respect  to  the 
d  n 

primary  measurement  n 


Using  the  actual  data  set  and  its  uncertainties,  Equation  I-Bl  can  be  solved  for 
the  uncertainty  of  each  of  the  thermal  properties  (Table  I-Bl). 
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TABLE  I-Bl  SUMMARY  OF  RESULTS  OF  UNCERTAINTY  ANALYSIS 


MEASUREMENT 

NOMINAL  VALUE 

UNCERTAINTY 

L 

0.0257  m 

±0. 0005  m 

Q 

4.596*10®  J 

±0.025*10®  J 

YTi 

59. 44°C 

±0. 56°C 

P 

2.453  Mg/m3 

±0.016  Mg/m3 

/  (Tt-y 

10. 02°C-hr 

+0.56°C-hr 

PROPERTY 

PERCENT  UNCERTAINTY 

a 

4.5 

k 

5.9 

c 

5.6 
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APPENDIX  1 1 - A 


JUSTIFICATION  FOR  MODELING  A  THREE-LAYER  ASPHALTIC 
CONCRETE  PAVEMENT  AS  A  HOMOGENEOUS  BODY 


OBJECTIVE 

The  objective  of  this  sub-study  is  to  justify  the  use  of  a  homogeneous 
model  of  a  three- layer  asphaltic  concrete  pavement  by  demonstrating  that  the 
errors  incurred  in  the  approximation  are  negligibly  small.  This  objective  is 
accomplished  by  comparing  a  characteristic  temperature  history  given  by  a  simpli¬ 
fied  analytical  model  of  the  transient  heat  transfer  in  the  actual  composite 
pavement  with  the  analogous  temperature  history  given  by  a  homogeneous  model  of 
the  pavement  and  showing  that  these  results  are  virtually  identical. 

ANALYSIS 

In  order  to  facilitate  the  comparison  described  above,  the  pavement  is  as¬ 
sumed  to  be  a  semi -infinite  body  initially  at  a  uniform  temperature,  T.. .  At 
time  t  =  0,  the  surface  of  the  pavement  is  assumed  to  undergo  a  step  change  in 
temperature  from  T.  to  Tg.  (The  surface  temperature  is  maintained  at  T£  for  all 
t  >  0. )  The  transient  conduction  heat  transfer  is  then  computed  using  two  dif¬ 
ferent  models  of  the  pavement.  First,  the  pavement  is  modeled  as  a  composite 
body  consisting  of  two  different  homogeneous,  constant-property  layers  (repre¬ 
senting  the  wearing  and  binder  courses).  The  neglect  of  the  base  course  and  the 
use  of  a  semi-infinite  geometry  are  reasonable  assumptions  for  the  purposes  of 
this  side  study  because  of  the  extremely  low  thermal  diffusivity  of  asphaltic 
concrete  and  the  relatively  short  heating  times  employed  in  a  typical  recycling 
process.  Thus,  the  wearing  course-binder  course  interface  temperature  history 
computed  with  the  two-layer  model  is  assumed  to  be  representative  of  the  actual 
interfacial  temperature  history  in  an  asphaltic  pavement  during  the  heating 
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phase  of  a  recycling  process.  Next,  the  pavement  is  modeled  as  a  homogeneous, 
constant-property  body.  The  temperature  history  at  a  distance  from  the  surface 
in  the  homogeneous  body  equal  to  the  interface  distance  in  the  two- layer  body 
is  then  compared  with  the  interface  temperature  history  computed  with  the  two- 
layer  model.  The  deviations  of  the  former  temperature  history  from  that  in  the 
two-layer  body  are  assumed  to  be  indicative  of  the  errors  incurred  in  modeling 
an  actual  composite  pavement  as  a  homogeneous  body. 

As  outlined  above,  we  first  assume  the  pavement  is  a  semi-infinite  body 
occupying  the  half-space  x  >  -  L,  as  shown  in  Figure  II-A1.  The  body  is  com¬ 
posed  of  two  layers  of  different  properties.  In  the  region  (-  L  <  x  <  0),  the 
thermal  conductivity,  specific  heat,  density,  and  thermal  diffusivity  are  given 
by  k,  Cp  and  a^,  respectively.  The  temperature  in  this  region  is  denoted  by 
T^  and  the  properties  (k^,  c^,  p^,  a^)  are  assumed  to  be  constants.  Similarly, 
in  the  region  (0  <  x  <  »)  the  properties  are  given  by  k^,  c 2,  p 2  and  and  the 
temperature  is  denoted  by  It  is  also  assumed  that  there  is  no  contact  re¬ 
sistance  between  the  two  layers  in  the  pavement.  Carslaw  and  Jaeger  (9)  show 
that  the  time-varying  temperature  distribution  in  material  2  (the  binder  course) 
is  given  by 


T2(x,t) 

T 

s 


where: 


— — —  I 

a  +  1  n=0 


ERFC 


(2n+l)  +  V 


<>n 


(II-A1) 


2VFo 


„  =  fUEL 

kipici 
a.t 

-  =  Fourier  number  (dimensionless  time) 

L2 
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(II-A2) 


Fo  = 


Evaluating  Equation  (II-A1)  at  x  =  0  yields  the  dimensionless  temperature  history 
at  the  interface  between  the  wearing  and  binder  courses  as 


■Tj  (t)l 


T2(o,t) 


o-l  n 


ERFC 


a  +  1  n-0  o+l 


2n  +  1 
2  V  Fo 


(II -A3 ) 


2  layer  s 

For  a  homogeneous  body  material  1  is  identical  to  material  2,  o  =  1,  and 
equation  (A3)  reduces  to: 


Tt  (t) 

— -  =  ERFC 

T  HOMOG 


1 

2  V  Fo 


(II-A4) 


This,  of  course,  is  just  the  classical  result  for  the  temperature  history  at 
a  distance  L  from  the  surface  of  a  semi-infinite  homogeneous  body  that  has  under¬ 
gone  a  step  change  in  surface  temperature. 

Finally,  we  assume  that 

%  error  in  modeling  the\  /%  error  in  the  interface 

composite  pavement  as  a)  =  l  temperature  history 
homogeneous  body  / 


(i 

) 

TI  | 

(  )  : 

^  *  100% 

Ts 

2  LAYER 

Ts  HOMOG 

(II-A5) 

( 

?-> 

s 

2  LAYER  1 

1 

Equation  (II-A5)  was  evaluated  with  the  aid  of  a  digital  computer  using  Equa¬ 
tions  ( I I- A3 )  and  (II-A4).  The  results  are  given  in  Figure  II-A2  for  a  range  of 
Fourier  numbers  (dimensionless  heating  times)  that  correspond  to  the  range  of 
actual  heating  times  expected  in  recycling  processes.  The  results  are  displayed 
for  the  parameter  a  =  .9  and  1.1.  These  values  of  o  correspond  to  +  20%  differ¬ 
ences  in  the  property  group  (pCk)  between  the  wearing  and  binder  courses.  Such 
differences  are  quite  conservative  and,  therefore,  the  errors  encountered  in  an 
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actual  heating  process  should  be  smaller  than  the  limiting  values  computed  with 
Equation  (II-A5)  using  these  values  of  a.  Figure  II-A2  shows,  however,  that  the 
maximum  possible  error  is  only  5%,  with  a  time-average  error  much  closer  to  2% 
when  the  Fourier  number  ranges  to  its  maximum  possible  value.  In  light  of  the 
many  illdefined  parameters  that  affect  the  heat  transfer  during  a  typical  re¬ 
cycling  process,  such  errors  are  assumed  to  be  negligibly  small.  Thus,  we  may 
conclude  that  the  transient  heat  transfer  in  the  upper  regions  of  a  three-layer 
composite  pavement  may  be  adequately  portrayed  by  modeling  the  pavement  as  a 
homogeneous  body. 
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APPENDIX  II-B 


SOME  ALTERNATIVE  MODELS  FOR  THE  COMBINED  RADIATION-CONVECTION  HEAT 
TRANSFER  AT  THE  UPPER  SURFACE  OF  THE  PAVEMENT 

As  stated  earlier  in  this  report,  the  exact  relative  roles  played  by  con¬ 
vection  and  radiation  in  the  surface  heating  of  asphaltic  pavement  are  currently 
not  well-understood.  Investigators  such  as  Spall  (26)  have  claimed  some  degree 
of  success  in  modeling  the  surface  heating  of  asphaltic  pavement  as  a  purely 
radiative  process,  however,  such  minimal  evidence  is  far  from  conclusive. 

Three  other  possible  models  of  the  surface  heating  process  as  listed  below: 

Model  1  -  Radiation  and  Forced  Convective  Heating 

Heat  transfer  to  the  pavement  surface  is  by  radiation  from  the  heater  sur¬ 
face  and  by  forced  convection  from  the  hot  combustion  gases  (which  are  assumed 
to  be  in  direct  physical  contact  with  the  pavement). 

Model  2  ~  Radiation  and  Forced  Convective  Cooling 

Heat  transfer  to  the  pavement  is  by  radiation  from  the  heater  surface.  A 
forced  convective  flow  of  ambient  air  is  assumed  to  be  drawn  in  from  the  sides 
and  thereby  cool  the  pavement  surface. 

Model  3  Radiation  and  Combined  Forced  and  Natural  Convection 

Heat  transfer  to  the  surface  is  by  radiation  from  the  heater  surface.  Na¬ 
tural  and/or  combined  natural  and  forced  convective  flows  of  ambient  air  are 
assumed  to  cool  the  pavement. 
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APPENDIX  II-C 


DESCRIPTION  OF  THE  FINITE-DIFFERENCE  TECHNIQUE 


The  Euler  method  is  an  explicit,  technique  that  computes  the  temperature 
at  the  end  of  a  time  interval  At*  =  (t*)  ,,dnr,  th« 


-  (t*)  using  the 


derivative  evaluated  at  the  beginning  of  the  time  interval,  or 


(T*)  <v+1>  =  (T.)  (v) 


dT*  ,  . 

+  (v)  (At*) 

dt* 


Rewriting  this  equation  in  terms  of  the  nondimensional  variables  gives 


T  (v+1)  =  T(v)  +  (  __L  )  (v)  (At) 

1  1  dt 


This  method  was  chosen  because  of  the  extreme  simplicity  of  its  application. 


Although  it  is  well-known  that  explicit  methods  are  subject  to  unstable  numerical 
oscillations,  it  is  just  as  well-known  that  these  methods  are  completely  stable 
if  a  small  enough  time  interval  is  employed  in  the  calculations.  The  usual  ob¬ 
jection  to  explicit  methods  is  that  the  small  time  intervals  required  to  insure 
their  stability  increase  the  computational  times  with  respect  to  implicit  meth¬ 
ods.  This  shortcoming  of  explicit  methods  is  far  outweighed  in  the  present  in¬ 
vestigation  by  the  ease  of  application  of  explicit  techniques  since  the  total 
computaional  time  required  will  be  quite  small,  independent  of  the  numerical 
technique. 
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APPENDIX  II-D 
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